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ABSTRACT 


Goals  of  this  research  have  been  to  identify  physical  processes  that  determine  the  dynamics 
of  the  marine  cloud  layers  and  to  quantify  roles  of  turbulence,  convection  and  thermal  radiation  that 
play  in  formation,  dissipation  and  stability  of  the  marine  cloud  layers.  And  immediate  objectives  of 
the  research  are  to  advance  turbulence  models,  use  efficient  numerical  schemes,  develop  computer 
simulation  programs,  simulate  the  marine  cloud  layers  and  compare  computer  results  with  published 
experimental  data  on  the  marine  cloud  layers  so  as  to  yield  insights  into  the  cloud’s  physical 

processes. 

For  these  objectives,  two  theoretical  models,  using  the  second-order-turbulence  closure  and 
the  large-eddy  simulation  (LES),  respectively,  have  been  developed.  In  addition,  efforts  have  been 
made  to  develop  a  hybrid  model  that  is  based  upon  a  framework  of  the  LES  model  but  uses 
turbulence  values  predicted  by  a  second-order-closure  model.  This  hybrid  model  preserves  details 
of  turbulence  but  eliminates  the  need  of  continuous  evaluation  of  the  multi  dimensional  turbulence 
equations.  The  model  will  improve  simulation  efficiency  and  conserve  computational  resources.  An 
evaluation  of  the  model  has  shown  excellent  agreement  between  simulation  results  with  relevant 
experimental  data  retrieved  from  reliable  web  site  resources. 

In  this  report,  a  progressive  development  of  the  second-order-closure,  LES  and  hybrid 
turbulence  models  for  simulation  of  the  marine  cloud  layers  is  described  and  a  comparison  of 
theoretical  results  with  experimental  data  is  presented  to  yield  better  insights  into  the  cloud’s  physical 

processes. 
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1. 


INTRODUCTION 


Goals  of  this  research  have  been  to  identify  physical  processes  that  determine  the  dynamics 
of  the  marine  cloud  layers  and  to  quantify  roles  of  turbulence,  convection  and  thermal  radiation  that 
play  in  formation,  dissipation  and  stability  of  the  marine  cloud  layers.  And  immediate  objectives 
of  the  research  are  to  advance  turbulence  models,  use  efficient  numerical  schemes,  develop  computer 
simulation  programs,  simulate  the  marine  cloud  layers  and  compare  theoretical  results  with 
published  experimental  data  on  the  marine  cloud  layers  so  as  to  yield  insights  into  the  cloud’s 
physical  processes. 

Studying  the  multi  dimensional  cloud  layers  requires  a  turbulence  model  that  permits  high- 
resolution  simulation  of  turbulence  of  different  scales  and  domain  sizes.  For  high  resolution, 
detailed  second-order-closure  model  of  turbulence  may  be  used,  and  for  calculation  efficiency,  the 
large-eddy-simulation  (LES)  model  may  be  developed.  Two  computer  simulation  programs,  using 
the  second-order-closure  model  (Chi,  1996  and  1998a)  and  the  LES  model  (Chi,  1998b), 
respectively,  have  been  developed  under  this  research  effort.  Additional  efforts  have  been  made  to 
develop  a  hybrid  model  that  is  based  upon  a  framework  of  the  LES  model  but  uses  turbulence  values 
predicted  by  a  second-order-closure  model  (Chi,  1999a).  In  addition,  confidence  in  the  simulation 
model  is  established  by  comparing  simulation  results  with  experimental  data  retrieved  from  reliable 
web  site  resources  (Chi,  1999a  and  1999b). 

In  this  report,  a  progressive  development  of  the  second-order-closure,  LES  and  hybrid 
turbulence  models  for  simulation  of  the  marine  cloud  layers  will  be  documented  first.  Numerical 
experiments  will  be  carried  out  to  examine  calculation  efficiencies  and  characteristics  of  the  marine 
cloud  layers.  Experience  in  tracking  relevant  experimental  data  from  reliable  web  site  resources  will 
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be  described.  To  increase  insight  into  the  marine  cloud  layers,  the  computer  simulation  results  are 
compared  with  experimental  data. 

2.  THEORETICAL  MODELS  FOR  THE  MARINE  CLOUD  LAYERS 

For  high  resolution,  a  detailed  second-order-closure  model  of  turbulence  has  been  developed. 
For  calculation  efficiency,  a  LES  model  has  been  developed.  Details  of  these  models  have  been 
reported  in  several  papers  by  the  author  (Chi,  1996,  1998a,  and  1998b).  For  ready  references, 
reprints  of  these  three  papers  are  appended  to  the  report.  In  addition,  a  hybrid  turbulence  model  that 
is  built  upon  a  framework  of  the  LES  model  but  uses  turbulence  values  predicted  by  a  second-order- 
closure  model  has  been  used  to  predict  stability  of  multi  dimensional  marine  cloud  layers.  This 
hybrid  model  preserves  details  of  turbulence  but  eliminates  the  need  of  evaluating  the  multi 
dimensional  turbulence  equations,  and  It  improves  prediction  efficiency  and  conserve  computational 
resources.  In  addition,  prediction  results  have  been  shown  in  excellent  agreement  with  relevant 
experimental  data  retrieved  from  reliable  web  site  resources.  Presented  below  in  Section  2.1  is  a 
second-order-closure  model  for  calculating  the  turbulence  values  and  in  Section  2.2  is  the 
framework  of  a  LES  model.  In  Section  2.3,  a  hybrid  model  that  is  based  upon  the  LES  framework 
but  uses  the  second-order-closure  turbulence  values  is  described  to  predict  stability  of  the  marine 
cloud  layers. 

2.1  A  Second-Order-Closure  Model  of  Turbulence 

When  these  assumptions  are  made:  (1)  vertical  coordinate  z  is  in  the  dominant  direction  of 
turbulent  mixing  of  atmospheric  air,  (2)  at  far  above  a  sea  surface,  the  geostrophic  balance  is 
maintained,  (3)  velocity  and  temperature  of  vapor  and  liquid  moisture  are  in  equilibrium,  and  (4) 
Boussinesq  approximations  are  used,  the  conservation  equations  for  momentum,  enthalpy  and  total 
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moisture  of  atmospheric  air  can  be  written  as: 


=f(v-V)-— -W— 
dt  8  dz  dz 


dt  8  dz  dz 


80  _  wd@  awe  i  3Fr 
dt  dz  dz  p  C  dz 
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where  the  upper-case  dependent  variables  represent  the  mean  fields  and  the  corresponding  lower¬ 
case  variables  represent  turbulent  fluctuations  of  the  same  variables.  Thus,  U,  V  and  W  are  mean 
velocity  components  in  the  x,  y  and  z  directions  and  at  the  time  t.  (Ug,Vg)  are  geostrophic  wind 
components  in  (x,y)  directions.  0  is  the  mean  moist-air  potential  temperature  which  is  defined  as 
[T-T0+(gz+LQv)/Cp] ,  and  Q  the  mean  total  moisture  mixing  ratio  which  is  defined  as  (Qv+Qc).  L  is 
the  water  latent  heat  of  vaporization,  Cp  the  constant  pressure  specific  heat,  T  the  temperature,  T0the 
referenced  temperature  at  z  equal  to  zero,  Qv  the  water  vapor  mixing  ratio,  Q  f  the  liquid  water  mixing 
ratio,  g  the  gravitational  acceleration,  f  the  Coriolis  parameter,  @v  the  virtual  dry  potential 
temperature  which  is  defined  as  [T(l+1 .609@v-Q)-To+gz/Cp],  q  the  standard  density,  v  the  kinematic 
viscosity,  and  P  the  buoyancy  coefficient.  The  second-order  correlations  in  these  equations  represent 
the  mean  turbulent  fluxes  of  momentum,  enthalpy  and  moisture  fluxes. 

Using  the  second-order-closure  assumption  (Mellor  and  Yamada,  1974  and  Moeng  and 
Arakawa,  1980),  transport  equations  for  calculating  the  turbulent  fluxes  may  be  written  as  follows: 
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In  equations  5  through  18,  the  time-derivative  terms  on  the  left-hand  side  model  the  transient 
variation  of  turbulence  correlations.  The  second-order  derivative  terms  on  the  right-hand  side  model 
turbulent  diffusion.  While  the  production  of  turbulence  due  to  buoyancy  is  modeled  by  terms  with 
buoyant  coefficient  P,  the  production  of  turbulence  due  to  friction  is  modeled  by  products  of  second- 
order  correlations  and  gradients  of  mean  variables.  Terms  with  coefficients  B,  C,  E  and  F  represent 
the  turbulent  redistribution.  The  characteristic  length  scale  X  is  equal  to  the  value  of  the  Blackadafs 
or  the  diffusion-length  scale  -  XB  or  XD-  whichever  is  the  smallest.  Turbulent  dissipation  is  modeled 
by  terms  with  coefficient  D.  Coefficients  k.  A,  B,  C,  D,  E  and  F  have  been  determined  semi- 
empirically  (Chi,  1994),  having  the  values  equal  to  0.35,  0.21,  0.46,  0.053,  0.132,  0.44,  and  0.23, 
respectively.  Numerical  procedures  and  computer  programs  have  been  developed  to  solve 
equations  1  through  18  with  appropriate  boundary  conditions  (Chi,  1996  and  1998a).  Prediction 
results  using  the  second-order-closure  model  will  be  presented  and  discussed  in  Section  4.1  of  this 

report. 
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2.2  A  LES  Model  for  Marine  Cloud  Layers 


When  these  assumptions  are  made:  (1)  the  Coriolis  force  is  negligible,  (2)  velocity  and 
temperature  of  vapor  and  liquid  moisture  are  in  equilibrium,  and  (3)  Boussinesq  approximations 
are  used,  the  conservation  equations  for  momentum,  enthalpy,  total  moisture  and  turbulent  energy 
of  atmospheric  air  can  be  written  as: 
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In  above  equations,  E  is  the  turbulent  energy  defined  as  q2/2,  and  the  buoyancy  terms  associated  with 
the  virtual  dry  potential  temperature  0V  have  been  defined  as  follows: 


^^=—+(0.609^-1)  —  —  For  Clear  Air  Layers  (25) 
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The  rate  of  dissipation  within  the  grid  volume  and  the  subgrid  eddy  coefficient  may  be  parameterized 
through 


e=0.\9Em/X 

(28) 

and  the  eddy  difiusivity  coefficients  have  been  defined  as  follows: 
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where  A  is  a  minimum  of  the  Blackadar’s  length  AB,  difiusion  length  AD, 
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And  the  turbulent  Prandtl  numbers  for  enthalpy  and  moisture  diffiision  (ag  and  au  )  are  both  equal  to 
0.75.  Numerical  procedures  and  computer  programs  have  been  developed  to  solve  equations  19 
through  31  with  appropriate  boundary  conditions  (Chi,  1998b).  Prediction  results  using  the  LES 

model  will  be  presented  and  discussed  Section  4.2  of  this  report. 

In  Section  2.3,  a  hybrid  model  that  retains  high  resolution  with  the  second-order-closure 
model  described  in  Section  2.1  and  calculation  efficiency  with  the  LES  model  described  in  Section 

2.2  will  now  be  presented. 

2.3  a  Hybrid  Turhnience  Model  for  Marine  Cloud  Lavers 

Second-  and  higher-order  turbulent  models  have  succeeded  in  advancing  theoretical 
understanding  of  the  marine  cloud  layers.  Complexity  of  those  models  often  makes  long-term 
simulation  of  the  cloud  layer  over  an  extensive  period  prohibitively  expensive.  A  hybrid  model  that 
retains  calculation  efficiency  of  the  LES  model  and  turbulence  details  of  the  second-order-closure 
model  has  been  developed.  The  hybrid  model  uses  Equations  19  through  23  from  the  LES  model  for 
conservations  of  momentum,  enthalpy,  moisture  and  mass,  respectively. 
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But,  instead  of  the  turbulent-kinetic-energy  Equation  24,  the  second-order-closure  model  equations 
1  through  18  will  be  used  in  evaluating  the  eddy-diflusivity-coefficient  K  values  defined  by  Equations 
29  through  32.  Numerical  Procedure  and  computer  algorithms  (described  in  Section  3  below)  have 
been  developed  to  evaluate  the  hybrid  model.  Prediction  results  using  the  hybrid  model  will  be 
presented  in  Section  4.3  and  verified  by  comparison  with  experimental  data  in  Section  5.2. 

3.  NTTMFRTCAI,  PROCEDURE  AND  COMPUTER  ALGORITHM, 

The  hybrid  model  presented  in  Section  2.3  uses  the  turbulence  equations  1  through  18,  the 
conservation  equations  19  through  23,  and  the  diffusivity  coefficients  defined  by  Equations  29 
through  32.  Numerical  procedures  and  computer  algorithms  have  been  developed  to  simulate  the 
marine  cloud  layers.  It  can  be  observed  above  that  conservation  and  turbulence  transport  equations 
1  through  23  described  above  can  be  written  in  the  following  general  form: 


da  +  yd  a  +  jyda  _  _  _d_^  da  ^  _  _d_(p  da )  +$ 
dt  dx  dz  dx  dx  dz  dz 
within  Domain  A 


(37) 


The  boundary  and  initial  conditions  may  be  written  as: 

—+aa+b=0,  On  boundary 
dn 


(38) 
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a(x,z,t=0)=ar](x,z).  At  initial  time  zero 


(39) 


In  above  equations,  a(x,z,t)  values  are  variable  values  of  velocity  components  U  and  W,  potential 
temperature  0,  and  moisture  mixing  ratio  Q;  <j>(x,z,t)  values  are  the  source/sink  strengths  for  a. 
Items  in  equation  37  describe  the  variables’  time  derivative,  convection,  diffusion  and  source/sink 
strength,  respectively,  'a'  and  V  values  in  equation  38  are  used  to  define  the  appropriate  boundary 
conditions. 

A  finite-volume  difference  scheme  (Patanka,  1980)  was  used  for  descritization  of  variables. 
The  flow  domain  to  be  simulated  was  divided  into  small  rectangular  elements;  shown  in  a  figure 
below  are  examples  of  several  such  elements.  It  can  be  seen  in  the  figure  that  variable  values  at  node 
points  of  the  elements  are  strategically  located.  The  node  point  for  velocity  value  is  located  in  the 
middle  of  the  rectangular  edge  and  the  node  point  for  other  variables  is  located  in  the  center  of  the 
rectangular  box. 

A  Figure  Showing  an  Element  with  Neighboring  Variables 
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Using  the  velocity  values  and  values  of  a  shown  in  the  figure  and  employing  an  upwind  logic, 
the  conservation  law  may  be  applied  to  obtain  an  expression  for  the  value  of  a  at  the  node  point  P 
in  terms  of  the  a-values  at  its  neighboring  node  points: 
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apap=aIfi.E+ay{p.w+a1a.T+aBaB+b 

(40) 

where  coefficient  values  ‘a’  and  ‘b’  may  be  expressed  in  terms 

of  the  difiiision  and  flux  parameters 

defined  as  follows: 
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In  addition,  in  solving  U  and  W  values  using  equation  40,  the  initial  pressure  values  will  have  to 
be  the  estimated  values;  consequently,  the  equation  40  will  yield  initially  approximate  and  Ww 

values.  Improved  U  and  W  values  may  be  calculated  from  their  approximate  values  using  equations: 

(47) 

U^Uw,app+—(P^'Pp) 

(48) 
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wr^^-pr) 

(49) 

(50) 

Substituting  U  and  W  values  in  equations  47  to  50  into  the  continuity  equation  23  yields  a  set  of 
linear  equations  for  pressure  values  at  the  solution  domain’s  node  points;  they  may  be  used  to  solve 
for  improved  pressure  values.  So  the  process  may  be  repeated  to  iterate  alternatively  for  improved 
values  of  velocity,  pressure  and  other  variables  at  the  node  points. 

Using  the  linear  equations  discussed  above  for  variable  values  U,  W,  ©,  Q,  E  and  P  at  node 
points,  a  computer  program  has  been  written  to  simulate  dynamics  of  the  marine  cloud  layers. 

4.  COMPUTER  SIMULATION  RESULTS 

4.1  Second-Order-Closure  Model  Simulation  Results 

The  second-order-closure  model  of  turbulence  described  in  Section  2. 1  has  been  used  to 
predict  transient  exchange  of  moisture  in  the  sea/air  interface,  mixing  of  moisture  in  the  marine 
atmosphere,  and  formation  of  the  cloud  layer.  Graphs  plotted  in  figure  1  were  reported  in  a  previous 
paper  written  by  the  author  (Chi,  1996);  they  were  used  as  initial  conditions  for  this  study. 

The  initial  conditions  shown  in  figure  1  were  obtained  on  the  assumptions  that  the  geostrophic 
wind  in  the  x-direction  was  at  lOm/s,  and  Coriolis  parameter  f  is  equal  to  1.0x1  O'4  s'1,  the  potential 
temperature  was  chosen  to  be  well-mixed  at  290K,  and  the  total  moisture  mixing  ratio  corresponding 
to  100  percent  relative  humidity  at  top  of  the  marine  planetary  boundary  layer  (MPBL).  Just  above 
the  MPBL  top  at  height  equal  to  1  kilometer,  there  was  an  inversion  layer  of 400-m  thick,  in  which 


12 


4.10  -0.10  4.00  0.00  0.06 


Fig.  1:  Initial  Velocity-Components  and  Potential-Temperature  Values 

temperature  was  increasing  at  0.2  K/m  and  the  moisture  mixing  ratio  was  decreasing  at  0.008 
g/kg.m,  respectively.  The  sea-surface  temperature  was  allowed  to  vary  diumally  within  the  range  of 
280  to  290F.  Conditions  shown  in  figure  2  were  for  the  instance  when  the  sea  surface  was  at  290  K. 

For  the  simulation  run  using  initial  conditions  shown  in  figure  1,  it  was  assumed  that  the  sea 
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surface  temperature  was  raised  abruptly  to  and  then  maintained  at  298  K,  and  the  air  total  moisture 
ratio  at  the  sea  surface  was  maintained  at  the  saturation  state.  After  simulation  runs  over  a  period 

j 

of  forty  hours,  a  large  body  of  data  was  generated.  Plotted  in  figures  2  are  predicted  contours  of 
temporal  mean  physical-property  values  of  air  in  the  marine  atmosphere;  plotted  in  figures  3  are 
predicted  contours  of  Reynolds-stress  and  turbulent-flux  values  in  the  marine  atmosphere.  Many  j 

interesting  characteristics  of  mixing  in  the  simulation  domain  can  be  observed  in  these  contours. 

Firstly,  rapid  interaction  at  the  air/sea  interface  can  be  observed  during  the  initial  period  of  zero  to 
600  minutes.  It  is  followed  by  a  calmer  development  of  the  marine  planetary  boundary  layer  for  about 
600  minutes.  Then,  during  the  next  600  minutes  transfer  of  enthalpy  and  moisture  continues,  as  can 
be  observed  from  the  predicted  contours  of  the  turbulent  thermal  and  moisture-flux  values  shown  in 
figures  2C  and  2D.  Also,  can  be  observed  in  Fig.  2D  is  the  formation  of  cloud  starting  at  the  200th 
minute,  rapid  deepening  of  the  cloud  layer  during  the  second  interval  of  200  to  1200  minutes,  and 
slower  growth  of  the  cloud  layer  during  the  period  of  1200  to  1800  minutes.  Finally,  a  steady  state 
is  established  at  around  the  2400th  minute. 

Plotted  in  figures  4  and  5  are  the  predicted  steady-state  conditions  at  the  end  of  the  simulation  i 

period  of  the  40th  hour.  The  graphs  shown  in  figure  4  can  be  compared  with  those  in  figure  1  for  the 
initial  conditions.  Changes  that  have  been  made  during  these  forty  hours  can  be  observed.  Warm  and 
moist  air  at  a  sea  surface  has  resulted  in  an  unstable  mixing  layer.  Variations  of  both  the  mean-  1 

quantity  and  the  turbulent-flux  profiles  can  be  observed  in  figure  4.  Specifically,  if  mean  horizontal 
velocity  components  plotted  in  figures  1 A  and  turbulent  kinetic  energy  values  plotted  in  figure  IB 
are  compared  with  those  corresponding  values  plotted  in  figure  4A  and  4H,  thickening  of  the 
unstable  boundary  layer  due  to  intense  turbulent  exchange  can  be  observed  in  figure  4A  and  4H. 
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Fig.  2:  Predicted  Contours  of  Mean  Atmospheric  Quantities 
2A  -  Mean  Wind  Velocity,  m/s,  2B  -  Mean  Potential  Temperature,  K, 

2C  -  Mean  total  moisture  mixing  ratio,  g/kg,  2D  -  Mean  liquid  moisture  mixing  ratio,  g/kg 


Fig.  3:  Predicted  Contours  of  Turbulent  Fluxes 
3A  -  Turbulent  kinetic  energy,  mA2/sA2,  3B  -  Principal  Reynolds  stress,  mA2/sA2 
3C  -  Turbulent  thermal  flux,  m.K/s,  3D  -  Turbulent  moisture  flux,  m.kg/kg.sA3 
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Fig.  4:  Predicted  Mean  Quant:  Ijs  and  Turbulent  Fluxes 
4A  -  Mean  velocity;  4B  -  Reynolds  stress; 

4C  -  Potential  temperature;  4D  -  Thermal  flux; 

4E  -  Total  moisture  ratio;  4F  -  Moisture  flux; 

4G  -  Liquid  moisture  ratio;  4H  -  Turbulent  energy. 
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Fig.  5:  Initial  steady  state  at  the  zeroth  hour 


Horizontal  in  kilometer* 


Fig.  6:  Contours  of  superposed  stream  functions  with 
entrainment  at  the  top 
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4.2  LES  Model  Simulation  Results 

The  LES  model  presented  in  Section  2.2  has  been  used  to  study  stability  of  the  marine  cloud 
layers.  Graphs  plotted  in  figure  4  and  5  predicted  by  a  second-order-closure  model  presented  in 
Section  4.1  above  can  be  used  as  a  set  of  initial  values  in  this  numerical  experiment.  It  may  be  noted 
that  solutions  for  this  set  of  graphs  were  obtained  by  using  the  main  stream  wind  velocity  U  equal  to 
10  m/s.  To  investigate  effects  of  cloud-top-warm-air  entrainment  on  stability  of  the  cloud  layer 
shown  in  figure  5,  a  stream  function  shown  in  figure  6  will  be  superimposed  to  the  main  flow.  In 
generating  the  stream  function  shown  in  figure  6,  the  top-down  entrainment  was  assumed  to  be  at  a 
maximum  rate  of  6  cm/s  is  at  the  top-left  comer  and  the  rate  was  reduced  sinusoidally  to  zero  at  the 
top-right  comer.  In  addition,  it  is  assumed  that  potential  temperature  of  the  entrained  air  is  at  five 
degrees  centigrade  higher  than  that  of  the  cloud  at  the  top. 

Using  the  initial  conditions  and  entrainment  rates  described  above,  a  LED  simulation  computer 
program  has  been  run.  Shown  in  figure  5  are  contours  of  the  initial  steady-state  potential  temperature 
values,  total  moisture  mixing  ratio  values  and  liquid  moisture  Mixing  values,  respectively.  Predicted 
dynamic  responses  of  the  cloud  layer’s  liquid  moisture  content  to  the  warm-air  entrainment  at  the 
top  are  shown  in  figure  7.  From  snapshots  shown  in  this  figure  of  the  cloud  contours  at  different 
times  (i.e.,  at  half,  one,  five  and  ten  hours  from  the  start  of  the  simulation  run),  dissipation  of  the 
cloud  layer  can  be  observed. 

A  comparison  of  this  example  with  that  presented  in  Section  4.1,  it  can  be  observed  that  the 
LES  model  is  superior  to  the  second-order-closure  model  in  calculation  efficiency.  However,  details 
of  turbulence  are  lost  in  trading  for  this  efficiency. 
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4.3  Hybrid  Model  SinmiatfoiLBesuHs 


In  order  to  preserve  details  of  the  second-order-closure  turbulence  and  calculation  efficiencies 
of  large-eddy  simulation,  a  computer  program  using  the  hybrid  model  developed  in  section  2.3  has 
been  written  to  study  stability  of  the  marine  cloud  layer.  Numerous  runs  of  the  computer  program 
have  been  made  to  simulate  dynamic  responses  of  marine  cloud  layers  under  a  variety  of  conditions. 
To  facilitate  later  comparison  of  simulation  results  with  TOGA-COARE  data  (Tropical  Ocean  Global 
Atmosphere  -  Coupled  Ocean  Atmosphere  Response  Experiment)  documented  by  Tao  and  Simpson 
(1993),  Miller  and  Riddle  (1994)  and  Lin  and  Johnson  (1996),  computer  runs  have  been  made  with 
internet  retrieved  experimental  stream-function  and  boundary-condition  values  as  input  data.  As  an 
example,  figure  8  shows  an  initial  experimental  dataset  for  the  stream-function  values,  and  figure  9 
shows  snapshots  of  cloud  profiles  predicted  by  the  hybrid  model,  at  several  different  instances, 
starting  from  the  zeroth  hour  to  the  24th  hour. 


Fig.  8:  Contours  of  a  set  of  retrieved  experimental  stream  function  values 


2  3 
DISTANCE  IN  lOOKM 


Fig.  9:  Snapshots  of  cloud  contours  predicted  by  a  hybrid  model  of  turbulence 
(Starting  from  an  initial  experimental  dataset  at  the  zeroth  hour) 
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5.  DATA  TRACKING  AND  THEORETICAL  VERIFICATION 
5.1  Internet  Atmospheric  Data  Tracking 


Several  data  sources  such  as  Goddard  Cumulus  Ensemble  (GCE),  GDACC  (Goddard 
Distributive  Active  Archive  Center),  FIFI  (First  ISLCP  Field  Experiment)  and  TOGA-COARE 
(Tropical  Ocean  Global  Atmosphere  -  Coupled  Atmosphere  Response  Experiment)  were  reviewed 
to  determine  the  time  interval,  grid  size  and  altitude  criteria.  It  was  determined  that  we  should  secure 
data  with  time  and  grid  intervals  at  six  hours  and  one  mile,  respectively.  Two  sources  of  data  TOGA- 
COARE  (Webster  and  Lukas,  1992)  and  GCE  (Tao,  1993)  were  then  considered. 

TOGA-COARE  DATA:  The  scientific  goals  of  COARE  are  to  describe  and  understand:  (1) 
the  principal  processes  responsible  for  the  coupling  of  the  ocean  and  the  atmosphere  in  the  western 
Pacific  warm-pool  systems,  (2)  the  principal  atmospheric  processes  that  organize  convection  in  the 
warm-pool  region,  (3)  the  oceanic  response  to  combined  buoyancy  and  wind-stress  forcing  in  the 
western  Pacific  warm-  pool  region,  and  (4)  the  multiple-scale  interactions  that  extend  the  oceanic  and 
atmospheric  influence  of  the  western  Pacific  warm-pool  system  to  other  regions  and  vice  versa.  To 
carry  out  the  goals  of  TOGA  COARE,  three  components  of  a  major  field  experiment  have  been 
defined:  interface,  atmospheric,  and  oceanographic.  The  experimental  design  calls  for  a  complex  set 
of  oceanographic  and  meteorological  observations  from  a  variety  of  platforms  that  carry  out  remote 
and  in  situ  measurements.  The  resulting  high-quality  dataset  is  required  for  the  calculation  of  the 
interfacial  fluxes  of  heat,  momentum  and  moisture,  and  to  provide  ground  thruth  for  a  wide  range  of 
remotely  sensed  variables  for  the  calibration  of  satellite-derived  algorithms.  The  ultimate  objective 
of  the  COARE  dataset  is  to  improve  air-sea  interaction  and  boundary-layer  parameterizations  in 
models  of  the  ocean  and  the  atmosphere,  and  to  validate  coupled  models. 
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Internet  web-site  data  was  used  to  review  the  data  sets  of  TOGA-COARE.  Several  sites  for 

data  were  reviewed.  The  web  site,  http://kiwi  atmos.colostate.edu/scm/toga-coare.html  was  found 

to  contain  mean  data  over  the  TOGA-COARE  IFA  region,  which  were  thought  appropriate  for 
testing  the  model  Using  FTP  commands  the  data  files  were  transfetred  to  a  main  flame  computer. 
Using  uncompress  command  the  compressed  files  were  uncompressed  and  converted  to  ASCII  files 
and  — — to  a  UDC  workstation  via  electronic  mails.  However,  owing  to  their  coarse  grid  sizes, 

higher  resolution  datasets  are  required  for  the  model  refinement. 

fiTF  MODFT  S  DATA:  The  Goddard  Cumulus  Ensemble  Model  (GCE)  is  maintained  by 

Mesoscale  Atmospheric  Processes  Branch  (MAPB)  at  Goddard  Space  Flight  Center  (NASA/GSFC). 
Scientists  in  NASA/GSFC  continuously  maintain  and  upgrade  the  system  (Tao  and  Simpson,  1993), 
and  the  GCE  model  is  simulated  and  operated  by  GSFS  system  analysts  and  operators  (Private 
communications).  The  GEC  model  may  be  used  to  interpolate  experimental  data  with  great  degrees 
of  accuracy.  The  following  procedures  were  followed  to  obtain  the  data  from  both  sources. 

Super  computer  CRAY  at  NASA  was  used  for  obtaining,  organizing,  processing  the  data. 
The  data  files  obtained  in  this  process  were  transmitted  to  the  UDC  workstation  via  electronic  mail. 
For  GCE  Model  data,  scientists  at  GSFC  provided  the  author  and  his  co-workers  at  UDC  with  the 
database  and  a  Fortran  program  including  a  list  of  parameters  via  electronic  mail.  Fortran  programs 
were  further  developed  and  modified  the  data  into  appropriate  formats  for  the  files.  The  data  from 
GCE  model  were  processed  using  the  Fortran  programs  in  UNIX  environment.  Then.  They  were 

transmitted  to  the  UDC  workstation  in  two  stages. 

Thp  First-stage  data  transmitted  included  the  following  parameters: 

parameter  (nx=512,nz=22) 
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parameter  (nhr=48)  !  number  of  hours  model  simulations  (h) 

parameter  (interval=2 1600)  !  Interval  to  read  in  data  (s) 
parameter  (ntimes=nhr*3 600/interval) 


real  latent  (nx,ntimes) 
real  pO  (nz) 
real  qcg(nx,nz,ntimes) 
real  qci(nx,nz,ntimes) 
real  qcl(nx,nz,ntimes) 
real  qcs(nx,nz,ntimes) 
real  qm(nx,nz,ntimes) 
real  rad_lw(nx,  nz,  ntimes) 
real  rad _sw(nx,  nz,  ntimes) 
real  sensible(nx,ntimes) 
real  temp(nx,nz,ntimes) 
real  u(nx,nz,ntimes) 
real  v(nx,nz, ntimes) 
real  vap(nx,nz,ntimes) 
real  w(nx,nz,ntimes) 
realzl(nz) 
real  z2(nz) 


!  surface  latent  heat  flux  (W/m**2) 

!  air  pressure  (mb) 

!  graupel  mixing  ratio  (g/kg) 

!  cloud  ice  mixing  ratio  (g/kg) 

!  cloud  water  mixing  ratio  (g/kg) 

!  snow  mixing  ratio  (g/kg) 

!  rain  water  mixing  ratio  (g/kg) 

!  long  wave  radiative  heating  rate  (K/hr) 

!  short  wave  radiative  heating  rate  (K/hr) 

!  surface  sensible  heat  flux  (W/m**2) 

!  air  temperature  (K) 

!  u-wind  speed  (m/s) 

!  v-wind  speed  (m/s) 

!  water  vapor  mixing  ratio  (g/kg) 

!  w-wind  speed  (m/s) 

!  model  grid  height  at  zl  levels  -  for  most  variables  here  (m) 
!  model  grid  height  at  z2  levels  -  where  2  is  (m) 


nx=512  (512  miles,  data  grid  =  one  mile) 
nz=22  ->  pressure  levels 
ntimes=32  (6-hrs  intervals) 


The  Second-stage  data  transmitted  included  the  following  parameter: 


real  tke(nx,nz) 

real  du(nx,nz) 

real  dv(nx,nz) 

real  dw(nx,nz) 

nz=22  ->  pressure  levels 

ntimes=32  (6-hrs  intervals) 


!  turbulent  kinetic  energy  (cm*cm/s/s) 

!  u-momentum  turbulence  rate  (cm/s/s) 

!  v-momentum  turbulence  rate  (cm/s/s) 

!  w-momentum  turbulence  rate  (cm/s/s) 


The  database  obtained  from  the  GCE  model  was  large.  Eight  six-hours  data  were  organized  in  eight 
different  data  files,  each  file  containing  one  six-hour  interval  of  data.  The  data  files  were  transmitted 


via  electronic  mails. 
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TtATA  PRESENTATION:  For  graphic  presentation,  they  were  further  decoded  at  UDC 
by  several  Fortran  programs,  and  were  translated  into  formats  recognizable  by  data  plotting  software 
packages  -  Surfer  and  Grapher  -  which  are  available  in  the  author’s  laboratory.  Presented  below  are 
plotted  sample  data  values  over  a  domain  with  a  height  of  3  -km  and  a  horizontal  distance  of  5 00-km, 

and  over  a  time  period  of  24  hours  at  12-hours  intervals: 

.  Contours  of  empirical  wind  flow  stream  functions  values 

.  Contours  of  empirical  total  moisture  mixing  ratio  values 

•  Contours  of  empirical  liquid  moisture  mixing  ratio  values 

•  Sea/air  interface  sensible  and  latent  heat  flux  values 


Fig.  10:  Contours  of  experimental  wind  stream-function  values 
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Fig.  1 1 :  Contours  of  experimental  total  moisture  mixing  ratio 


Fig.  12:  Contours  of  empirical  cloud  profiles 
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Fig.  13:  Sea  surface  sensible  and  latent  heat  flux  values 
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The  hybrid  model  developed  in  the  program  preserves  details  of  turbulence  of  the  second- 
order-closure  model  and  calculation  efficiency  of  the  large-eddy  simulation.  To  establish  confidence 
in  prediction  accuracy,  the  cloud  profiles  predicted  in  Section  4.3  using  the  hybrid  model  shown  in 
figure  14A  are  compared  with  the  empirical  dataset  presented  in  Section  5.1  plotted  in  figure  14B. 


DISTANCE  IN  100KM  DISTANCE  IN  lOOKM 


Fig.  14 A:  Predicted  Cloud  Profiles 


Fig.  14B:  Empirical  Cloud  Profiles 


6.  CONCLUSION 


For  goals  of  research  at  UDC  to  identify  physical  processes  that  determine  the  dynamics  of 
marine  cloud  layers  and  to  quantify  roles  of  turbulence,  convection  and  thermal  radiation  that  play 
in  formation,  dissipation  and  stability  of  the  marine  cloud  layers,  the  author  and  his  cc-workers  at 

UDC  have  endeavored  to  achieve  the  objectives: 

•  To  advance  turbulence  models  using  efficient  numerical  schemes, 

•  To  make  available  computer  simulation  programs  for  predicting  stability 
of  marine  cloud  layers, 

•  To  retrieve  web-site  data  of  marine  cloud  layers  from  reliable  sources, 

•  To  establish  confidence  in  models  through  experimental  verification,  and 

•  To  yield  insights  into  the  cloud’s  physical  processes. 

For  these  objectives,  two  turbulent  models,  using  the  second-order-closure  and  the  LES 
model,  respectively,  have  been  formulated.  A  finite-volume  procedure  has  been  developed  to  solve 
the  resultant  equations.  In  addition,  to  retain  turbulence  details  of  the  second-order  turbulence  and 
calculation  efficiencies  of  the  LES  model,  a  hybrid  model  has  been  assembled.  It  uses  a  multi¬ 
dimensional  framework  of  the  LES  model  but  uses  a  set  of  unidirectional  second-order  equations  to 
calculate  the  turbulent  intensity  values. 

For  experimental  verification  of  the  hybrid  model,  theoretical  results  have  been  compared  with 
relevant  empirical  data  retrieved  from  reliable  sources.  Several  data  sources  such  as  TOGA-COARE 
(Tropical  Ocean  Global  Atmosphere  -  Coupled  Ocean  Atmosphere  Response  Experiment),  GCE 
(Goddard  Cumulus  Ensemble),  FIFE  (The  First  ISLCP  Field  Experiment),  GDAAC  (Goddard 
DAAC)  and  others  were  evaluated  based  on  the  time-interval,  grid-size  and  domain  criteria. 


29 


According  to  these  criteria,  two  sources  of  data  TOGA-COARE  (experimental)  and  GCE  model 
(semi-empirical)  were  examined  in  detail.  Through  cooperation  with  the  NASA  scientists  and  a 
permission  of  using  their  CRAY  computer,  several  sets  of  the  GCE  data  were  decoded  and 
transmitted  to  the  author.  These  data  have  been  used  to  verify  the  UDC’s  model  simulation  results. 
Figure  14A  shows  snapshots  of  cloud  profiles  that  have  been  predicted  by  the  hybrid  model  under  the 
same  wind-stream-function,  air/sea  interface  sensible  and  latent  heating  fluxes  and  other  boundary- 
condition  values  that  were  provided  by  the  empirical  data.  Although  predictions  are  not  expected  to 
have  a  complete  agreement  with  the  experimental  data  (as  precipitation  is  not  considered  in  the 
present  theory),  a  good  degree  of  similarity  can  be  observed  between  the  predicted  cloud  profiles 
shown  in  figure  14A  and  those  empirical  cloud  profiles  shown  in  figure  14B. 
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ABSTRACT 

A  study  has  been  made  to  identify  physical  processes  that 
determine  the  dynamics  of  the  marine  atmosphere  and  its  cloud  layers. 
The  effort  starts  from  formulation  of  die  governing  equations  for 
conservation  of  momentum,  enthalpy  and  moisture  of  atmospheric  air. 
The  turbulence  transport  equations  are  derived  to  calculate  the 
Reynolds-stress  and  turbulence-flux  correlations  that  appear  in  the 
conservation  equations.  The  thermal  radiation  equations  are  derived 
to  calculate  the  radiation  flux  divergence.  A  virtual  dry  potential 
temperature  value  is  introduced  to  account  for  the  effects  of  water 
vapor  condensation  and  cloud  formation  on  the  transport  of  the 
turbulence  quantities.  A  simulation  computer  program  is  developed  to 
calculate  dynamic  responses  of  the  marine  atmosphere  and  its  cloud 
layers.  Results  of  numerical  experiments  have  led  us  to  a  better 
understanding  of  die  effects  of  heat  and  moisture  exchange  between 
air  and  sea,  cloud-top  radiative  cooling,  and  interactions  of 
atmospheric  turbulence  and  thermal  radiation  on  the  dynamics  of 
marine  cloud  layers. 


CONSERVATION  AND  TURBULENCE  EQUATIONS 
Let  U,  V  and  W  be  the  mean  velocity  components  in  the  east,  north 
and  vertical  (x,y,z)  directions  and  at  the  time  t,  FR  the  radiation  dux, 
0  the  mean  moist-air  potential  temperature  defined  as 
[T-Kgz+LQ.yCpJ,  and  Q  the  mean  total  moisture  mixing  ratio  of  vapor 
and  liquid  water  (Qv-Ki).  The  lowercase  symbols  u,  v,  w,  0  and  u> 
represent  the  turbulent  fluctuation  quantities  corresponding  to  their 
respective  mean  quantities  U,  V,  W,  ©  and  Q.  In  addition,  the 
following  definitions  will  be  used:  L  for  the  water  latent  heat  of 
vaporization,  Cp  the  constant  pressure  specific  heat,  T  the  temperature, 
g  the  gravitational  acceleration,  f  the  Coriolis  parameter,  0,  the  virtual 
dry  potential  temperature  defined  as  [T(l+i.609QT-Q)+gz/Cp),andp 
the  density.  Far  above  a  sea  surfkcc,  the  geostrophic  balance  is 
assumed,  h  gives  the  following  relationship  between  the  pressure  field 
and  horizontal  velocity  components  Ug  and  Vg  (Holton,  1972): 


1  dP 


par 


(1) 


INTRODUCTION  j  dp 

Advances  in  turbulence  modeling  and  computational  techniques  — —  (2) 

have  made  it  a  common  practice  to  study  numerically  die  turbulent-  P  hV 

flow  problems  (Mdlor  and  Yamada,  1974;  Moeng  and  Aiakawa, 

19S0;  Moeng  and  Randall,  1984;  Chi,  1991).  Several  studies  of 

turbulent  flows,  tg.,  using  the  eddy-viscosity  model  for  a  tomado-like  h  addition,  equilibrium  of  velocity  and  temperature  between  phases 

vortex  (Chi  and  Jih.  1974),  the  large-eddy-simulation  model  for  end-  is  assumed,  and  a  Boussinesq  approximation  is  used.  Consequently, 

wall  boundary  lawn;  of  intense  vortices  (Chi,  1977),  the  e-k  model  for  *e  conservation  equations  for  momentum,  enthalpy  and  moisture  of 

vortex  flows  over  the  water  surface  (Chi,  1987)  and  the  second-order  atmospheric  air  can  be  written  as: 

closure  model  for  heat  and  moisture  transfer  in  marine  atmospheres 
(Chi,  1994),  have  been  reported  earlier.  In  this  paper,  a  study  of  the 
simultaneous  effects  of  turbulent  heat  and  moisture  transfer,  water 
vapor  condensation  and  thermal  radiation  on  the  dynamics  of  marine 
cloud  layers  is  presented. 
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Here,  Ej  represents  the  emissi vity  value  E  at  the  height  z  where  the  FR 
value  is  being  calculated. 

(4) 

NUMERICAL  PROCEDURE  AND  COMPUTER 
ALGORITHM 

The  conservation  and  turbulence  transport  equations  discussed 
above  can  be  written  in  the  following  general  form: 

(5) 

§!!L+U—--^-/y—)+<i>=0,  within  domain  Q  (10) 
dt  JdXj  &/  dx/  V 

(6) 

and  their  boundary  and  initial  conditions  may  be  written  as: 


The  turbulent  transport  equations  (Chi,  1991)  derived  from  the 
second-order-closure  assumption  can  be  used  to  solve  a  set  of  the 
fifteen  turbulent  variables  u7,  v7,  w7,  TTv,  uw,  v  w,  u"ff,  T6,  wff, 
uw,  vw,  wu>,  TP,  and  w7;  they  include  the  variables  wu,  wv, 
wS*,  and  wuin  Eq.  (1-4).  To  allow  for  water  vapor  condensation  and 
doud  formation,  the  buoyancy  terms  associated  with  fluctuation  of  the 
virtual  dry  potential  temperature,  vUy,  wff  „and  wF  ¥ 
appear  in  file  turbulence  transport  equations;  they  may  be  calculated 
from  the  predicted  turbulence  variables,  using  the  following  equations 
for  the  clear  and  cloudy  layers,  respectively: 

|5>S&+(0.609<M)t^  (7) 


gg>l*l-6P9^Ig-4>-ifo  (8) 
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Here  5  represents  any  turbulence  quantities,  u,  v,  w,  etc.;  4>  is  equal  to 
CpT/L;  y  ts  equal  to  (L/C  p)cQ7dT;  and  superscript  *  represents  values 
at  the  saturation  state. 


THERMAL  RADIATION  EQUATIONS 

The  radiation  calculation  requires  a  large  amount  of  time.  Simple 
emissi  vity  approach  using  the  gray-gas  assumption  has  been  used. 
The  theory  discussed  by  Sparrow  and  Cess  (1966)  was  followed.  By 
applying  the  boundary  conditions  of  the  blackbody  radiation  at  the  sea- 
surface  and  zero  downward  radiation  at  the  top  of  the  atmosphere,  and 
letting  E  be  the  mixed  water-vapor-and-cloud  emissivity  value,  the 
thermal  radiation  flux  FR  can  be  calculated  by  the  equation: 

Fr  =oT*e  f  zoT*e  ^^dz 
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y—+aa+b=0,  on  boundary  T  (11) 
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a(xyj,0 )=a0(xyj:),  at  initial  time  zero  (12) 


where  a(z,t)  are  mean  values  of  velocity  components,  enthalpy, 
moisture  and  turbulent  correlations;  4>(z,t)  values  are  sources/sinks  for 
a;  the  Four  items  in  Eq.  (10)  describe  the  system’s  inertia,  convection, 
diffusion  and  source/sink  strength,  respectively,  and  values  of ’a'  and 
V  in  Eq.  (1 1)  may  be  specified  to  enforce  the  appropriate  boundary 
conditions. 

A  finitedement  procedure  (Chi,  1994)  is  used  for  the  discretization 
of  spatial  variables.  The  flow  domain  to  be  simulated  is  divided  into 
small  finite  elements.  Multiplying  both  sides  of  Eq.  (10)  by  a 
weighing  function  o  and  integrating  over  the  finite  element  Q«,  the 
weighted  residual  equation  set  can  be  derived: 


f  a>[%+l7^-^-(Y^)+4>yt=0  (13) 

Jo.  dt  &. 


.  da  d  da 


Integrating  the  second-order  term  in  Eq.  (13)  by  parts  and  substituting 
Eq.  ( 1 1 )  for  normal  gradients  at  the  dements  boundary  result  in  the 
following  equation: 

r  [c,>j£  (aua+ix>)cfa=0  (14) 

i  o#  dt  dx.  cx.  1  dx.  -fr, 

The  weighing  functions  g>  and  dependent  variables  a(x,y,z,t)  are 
interpolated  in  each  dement  Q,  as  follows: 


co/z)=ilrn(z)Qn  (15) 

06) 
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where  T|f„  denotes  the  interpolation  functions  at  the  nth  node  of  each 
finite  dement,  and  A,  and  Q  „are  the  dement  nodal  values  of  the 
variable  c  and  the  weighing  function  u,  respectively. 

Substituting  Eq.  (15)  and  (16)  into  Eq.  (14)  yields  the  following 
finite-element  equation: 


Vi„=0 


where  C_  K_.  U„  end  Rw  »e=  mth-row  and  nth-duirn 
dements  of  square  matrices  which  describe  distributions  of  the  finite 
dements  capacitance,  difiusivity,  convectiveness  and  surface  flux, 
respectively,  and  L,  is  the  nth  dement  of  a  load  vector  for  the  finite- 
element  source/sink  terms.  ,  ,  . , 

Assembling  the  finite-element  Eq.  (17)  over  the  whole  flow  field 
y  !s  a  set  of  ordinary  differential  equations: 


+[K+ '  r+H]{A)  ♦(!]  =0 


(18) 


Where  the  vector  {A}  contains  all  node  point  variables  lying  within 
the low  domain  Q;  the  matrices  [q,  [K],  [U]  and  [H]  represent  global 
distnbutions  of  capacitance,  diffusiveness,  convection  and  surface 
conductance,  respectively,  and  [L]  is  the  global  load  vector. 

The  finite-element  procedure  described  above  has  trans-formed  the 
initial  value  partial  differential  conservation  and  turbulence  transport 
equations  into  a  large-order  system  of  the  ordinary  differential 
equation  set  (18).  Solution  of  (18)  may  start  from  a  Taylor  senes 
expansion  of  {A}  about  time  t 


ijxnj  % 

(A)  ,=tA  )+Af[6 - — *(l-8)— — — ]  (19) 

■  dt  dt 

where  0  is  an  implicit  coefficient  Multiplying  both  sides  of  Eq.  (19) 
by  capacitance  matrix  [C],  Eq.  (19)  may  be  re-wntten  as. 

(FlA..1))=(atlAl..1-lA.)-(A7)e^^-(A/)(l  (20) 

With  expressions  fortune  derivatives  provided  by  Eq.  (18),  Eq.  (20) 
becomes: 


iFI=[q..,(IA)..,-IA)J 

«♦! 

+Af(l  -e)([K>  ♦HIM  <21) 


The  equation  set  (21)  may  be  solved  for  {A}..,  values  at  time  V, 

from  the  known  {A},  values  at  timet,  by  a  Newton's  iterative  process. 

A  computer  program  has  been  written  to  facilitate  the  process.  The 
resultant  computer  program  can  be  used  to  simulate  the  dynamics  of 


the  marine  atmosphere  and  its  cloud  layers.  Presented  below  are  two 
examples  of  simulation  results. 


RESULTS  AND  DISCUSSION 

Firstly,  a  numerical  experiment  was  made  to  show  responses  of  a 
marine  atmosphere  to  variations  of  the  sea  surface  temperature.  Here 
the  geostrophic  wind  velocity  aligned  in  the  x  direction  issttati 
constant  value  of  10  m/s,  the  Coriolis  parameter  f  is  set  at  1-0x10  * »  1 
and  the  sea  surface  roughness  parameter  z,  is  set  at  0.1  m.  tatialty, 


mixed  at 290  *K  and  have  a  total  mixing  ratio  value  corresponding  to 
100  percent  relative  humidity  at  the  top  of  file  marine  planetary 
boundary  layer  (MPBL).  Just  above  the  MPBL  top  at  z  equal  to  1  km, 

there  is  an  inversion  layer  4(X)m  thick,  in  which  temperature  increases 

and  mixing  ratio  decreases  at  0.02  *K/m  and  SxlCT*  m'1,  i«P«tive>y. 
Above  the  inversion  top,  the  gradient  of  temperature  is  41.005  TC/m 
and  that  of  the  mixing  ratio  is  zero.  Subsequently,  the  sea-surface 

temperature  is  allowed  to  vary.  Numerical  simulations  are  made 
experiment  on  the  dynamic  responses  of  the  marine  atmosphere  to 
variations  of  the  sea-surface  temperature. 

Figure  1  shows  a  hypothetical  cyclic  sinusoidal  temperature 
variation  of  Ac  sea  surface.  Calculations  were  allo  wed  to  proceed  for 
five  days  with  Ac  sea-surface  temperature  repeating  cyclically  every 
twenty  four  hours.  Shown  in  Fig.  2  through  4  are  predicted  cyclic 
variations  over  a  twenty-four-hour  period  of  the  mean  wind-velocity 
and  air-potential-temperature  distributions  at  different  time.  Predicted 
responses  of  turbulent  energy,  horizontal  components  of  the  Reyndds 
ctr^c  and  vertical  comoonent  of  Ae  enthalpy  flux  are  shown  in  Fig. 


5  through  8.  .  , 

Many  interesting  characteristics  of  Ae  marine  planetary  boundary 
layer  (MPBL)  can  be  observed  in  Aese  figures.  It  can  be  seen  m  Fig. 

1  that  Ae  sea-surface  temperature  was  set  to  rise  during  Ae  hours  of 
six  to  eighteen  and  fall  during  Ae  hours  of  eighteen  to  six.  As  Ae  sea- 
surface  temperature  falls  (rises),  stability  of  Ae  boundary  layer 
increases  (decreases).  The  trends  can  be  seen  in  Fig.  2  and  3  by 
decreasing  in  Ae  marine  surface  boundary  layer  thickness  during  Ae 
hours  of  eighteen  to  six  and  increasing  in  Ae  marine  surface  boundary 
layer  thickness  during  Ae  hours  of  eighteen  to  six.  In  pmticular, 
instability  of  the  marine  thermal  boundary  layer  can  be  observed  in  Ae 
vicinity  of  Ae  cighteenA  hour  when  the  sea-surface  temperature  is  at 
its  maximum.  It  can  be  observed  also  in  Fig.  5  where  contours  of  Ae 
turbulent  energy  values  at  different  height  and  time  are  showm 
Instability  of  Ae  MPBL  in  Ae  vicinity  of  Ae  eighteen  A  hour  can  be 
seen  also  in  this  figure.  AlAough  it  can  be  observed  in  this  figure,  the 
maximum  thickness  of  Ae  MPBL  in  Ae  experiment  does  not  occur  at 
Ae  exact  time  of  Ae  maximum  sea-surface  temperature  at  the 
cighteenA  hour  (but  at  approximately  Ae  21st  hour)  because  of  delay 
in  Ac  responses.  Similarly,  Ae  effects  of  stability  on  distribution  of 
Ae  horizontal  components  of  Ae  Reynolds  stress  and  Ae  vertical 
component  of  the  heat  flux  can  be  seen  in  Fig.  6  to  8.  Agam, 
instabilities  of  Ae  MPBL  can  be  seen  in  Aese  figures  by  J* 

maximum  stress  and  flux  values  and  increasing  in  Ae  MPBL 
thickness  in  the  vicinity  of  Ae  eighteenA  hour. 

For  Ae  second  case  of  numerical  experiments,  its  initial  conditions 
at  the  zeroA  hour  are  obtained  from  the  eighteenth-hour  results  of  Ae 

first-case  experiments  presented  above,  In  this  case^he sea-surface 

temperature  is  changed  at  Ae  zeroA  hour  to  298*K,  ami  it  is 
maintained,  Aereafter,  at  this  constant  value.  The  ses-surface-air 
moisture  mixing  ratio  is  maintained  at  Ae  corresponding  saturated 
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Fig.  3  Y-Component  mean  wind  velocity,  V  m/s  Fig.  4  Mean  air  potential  temperature,  0  *K 
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Fig.  9  Contours  of  predicted  cloud  intensity  of 
the  liquid-water  mixing  ratio  in  g/kg 


state.  Plotted  in  Fig.  (9)  are  contours  of  the  predicted  cloud  intensity 
of  the  liquid-water  mixing  ratio  values  in  g/kg.  It  can  be  seen  in  this 
figure  that  the  cloud  starts  to  form  at  the  14th  hour.  Its  depth  grows 
firstly  downward  because  of  the  upward  transport  of  the  water  vapor 
from  the  sea  surface  and  the  upward  drop  in  the  static  temperature  of 
the  air.  From  the  25th  hour  onwards,  because  of  the  cloud-top 
radiative  cooling,  the  depth  of  the  cloud  increases  at  the  top.  Also,  it 
can  be  seen  in  this  figure  that  the  steady  state  of  the  cloud  layer  was 
reached  at  about  the  48th  hour. 


CONCLUSION 

In  conclusion  a  numerical  study  of  turbulence  in  marine 
atmospheres  has  been  made.  Results  of  predictions  have  led  us  to  a 
better  understanding  of  the  effects  of  air/sea  heat  and  moisture 
exchange,  cloud-top  radiative  cooling,  and  interactions  of  atmospheric 
turbulence  and  thermal  radiation  on  the  dynamics  of  marine  cloud 
layers. 
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^Goals  of  our  atmospheric  research  at  UDC  are  to  identity  physical 
processes  that  determine  the  dynamics  of  the  cloud  “J* 

quantify  the  roles  of  turbulence,  convection  and  thermal  radr  aho 
JSTpky  in  formation,  dissipation  and  stability  of  the  c^l  ^e 
Our  immediate  objectives  are  to  advance  theoreUcal  mode1^ 
efficient  numerical  schemes  and  develop  computer  pro^aMto 
simulate  the  marine  cloud  layers.  Comparison  of  computer  results 
with  published  field  observations  will  yield  insights 
layers’  physical  processes.  While  a  companion  paper  FEDSM98- 
4954  develops  a  large-eddy  turbulent  model  to  SUTlu'atefJb^J 
the  marine  cloud  layers,  this  paper  describes  a  second-ord« ’ 
turbulent  model  to  simulate  turbulent  mixing  of  the  moist 
atmospheric  air  in  general  and  the  marine  atmospheres  in  particular 

Of  hot,  moismr.  and 

environmentalists.  Advances  in  turbulent  modeling  “ 

better  understanding  of  climatic  phenonmena  A  -  cloud  stability 
tornado  motion,  and  severe-storm  formation  -  “d  “vuonmortal 
Qualities91013  -  air  pollution,  acid  deposition  and  global  warming. 
Slequently,  improvements  can  be  made  in  climatic  prediction  and 

environmental  control.  .  . 

Much  of  our  understanding  of  mixing  processes  in  atoospheres 
has  come  ftorn  careful  observation  and  sound  theoretical  modeling. 
We  it  is  yet  impossible  to  have  a  generalized  turbulace&eay 
JTunivemal  ph<!nomena,  semi-empirical  mode  s  w^ere* 

degrees  of  complexity  have  been  developed  with  confi^ncem 
shSSng  numerous  practical  phenomena:  the  mixing-length  theory 


has  been  used  to  model  turbulent  boundary  layers  on  flat  plates  '- 
_ J  model  has  been  used  for  simulating  the 


the  second-order  diffusion  model  has  oeen  usea  w  e  — 

SI*,  lsycrs,*  and 

io  study  the  meso-scale  turbulence  in  atmospheres.  Tins  andrf 
has  used  the  edAr-viseosity  model  to  simulate  a  tornado-like  vortex, 
SSSLbuleuce  model  for  end-wdl  boumtay  layers  of 
SX  die  e-k  model  for  vorm* 
surface  and  the  second-order  closure  model  for  the  manne  cloud 
layer.  It  is  intended  in  this  study  to  present  a  hybrid 
atmospheric  turbulence  that  uses  a  second-order  clo5^turb 
diffusion  model  for  simulating  atmosphenc  mixing  layers  and 
large-eddy  turbulent  model  for  simulating  meso-scale  entnunrnem 
inifoe  marine  atmosphere.  Presented  ini this  paper: is  - 
diffusion  model  for  simulating  moisture  nmong  m  fce mwnc 
atmosphere;  a  large-eddy  model  for  simulating  the  cloud-layer 
entrainment  will  be  presented  in  a  companion  paper. 

NOMENCLATURE 

C  =  specific  heat 

a!b  C  D  E  &F  values  in  turbulent  flux  equation 

’  =’  second-order  closure  equations'  empirical  constants 

f  =  geostrophic  Coriolis  coefficient 

• .  =  thermal  radiation  flux 

„R  =  gravitational  acceleration 

I  =  water-vapor  latent  heat  of  vaporization 

q  =  square-root  value  of  squared  mean  of  the  turbulent 

fluctuating  velocity 
t  =  time 

T  =  static  temperature  of  atmosphenc  air 

(u,v,w)  =  turbulent  fluctuating  velocities  in  (x,y,z)  ons 
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(U.V.W)  =mean  wind  velocity  components  in  (x,y,z)  directions 
(Ug,Vg)  =geostrophic  wind  components  in  (x,y)  directions 
(x,y,z)  =east-,  north-  and  vertical-direction  coordinates 
a  =a  generalized  variable  in  a  numerical  scheme 
P  =  buoyant  coefficient 

0  =  turbulent  fluctuating  value  of  0 

0V  =  turbulent  fluctuating  value  of  9V 
0  =  mean  potential  temperature  =  T-Kgz+LflJ/Cp 

eT  =  mean  virtual  potential  temperature 
=  T(l+l.G09fyQ}*&/C9 
k  =  Prandtl  mixing  length  constant 

A  =  turbulent  length  scale 

to  =  turbulent  fluctuating  value  of  fi 

Q  =  mean  total  moisture  ratio 

Q,  =  mean  liquid  moisture  ratio 

Q*  =  mean  vapor  moisture  ratio 

Overlined  values 

=  turbulent-flux  values 
Other  symbols  in  numerical  procedures 
=  those  defined  in  the  text 


(8) 

dt  dz  dz  dz  X 


(9) 

dt  dz  dz  oz  A 


dt  dz  dz  dz  dz  X 


dt  dz  dz  dz  dz  X 


dt  dz  dz  dz  A 


CONSERVATION  AND  TURBULENCE  EQUATIONS 

When  these  assumptions  are  made:  (1)  at  far  above  a  sea  surface, 
the  geostrophic  balance  is  maintained,  (2)  velocity,  temperature  of 
vapor  and  liquid  moisture  are  in  equilibrium,  and  (3)  Boussinesq 
approximation  is  used,  the  conservation  equations  for  momentum, 
enthalpy  and  total  moisture  of  atmospheric  air  can  be  written  as: 


dt  dz  dz  dz  A 


^±(AgX^)-™^-™^-EZ™ 
dt  dz  dz  dz  dz  X 


®L=KV-VV—-W— 

rit  ?  r)r 


2a=l(AgX&!t)-™&-™¥-El™ 

dt  dz  dz  dz  dz  X 


fSLjfjj  _U)Jn»_war 

dt  1  dz  dz 


dt  dz  dz  dz  A 


=  — 

dt  dz  dz  pC  dz 

p 


dtdz  dz  dz  A 


dQ  ^ydQ  _  dnti) 
dt  dz  dz 


dt  dz  dz  dz  dz  A 


In  above  equations,  values  for  overlined  turbulent  -flux  values  can 
be  calculated  by  the  turbulent  transport  equations  using  the 
assumptions  of  Mellor  and  Yamada’s  second-order  closure  model:7 

*»L,lAAgX*?L)-2™™  -B1(?-3L)-d£  (5) 

dt  dz  dz  dz  X  3  3A 


<6) 


^L,l(3AqX^.)^we,-B^w2-3l).D^. 
dt  dz  dz  A  3  3A 


In  above  equations,  A  is  a  characteristic  length  scale  that  is  equal  to 
the  value  of  the  Blackadar’s  or  the  diffusion-length  scale  -  AB  or  AD  - 
whichever  is  the  smallest 

In  equations  5  through  1 8,  the  time-derivative  terms  on  the  left- 
nand  side  model  the  transient  variation  of  turbulence  correlations. 
The  second-order  derivative  terms  on  the  right-hand  side  model 
turbulent  diffusion.  While  the  production  of  turbulence  due  to 
buoyancy  is  modeled  by  terms  with  buoyant  coefficient  p,  the 
production  of  turbulence  due  to  friction  is  modeled  by  products  of 
second-order  correlations  and  gradients  of  mean  variables.  Terms 
with  coefficients  B,  C,  E  and  F  represent  the  turbulent  redistribution. 
Turbulent  dissipation  is  modeled  by  terms  with  coefficient  D. 
Coefficients  k.  A,  B,  C,  D,  E  and  F  have  been  determined  semi- 
empirically,3  having  the  values  equal  to  0.35,  0.21,  0.46,  0.053, 
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0.132, 0.44,  and  0.23,  respectively. 

numerical  procedure  and  computer 
algorithm  ,  lfi 

The  conservation  and  turbulence  transport  equations  1  to  18 
described  above  can  be  written  in  the  Mowing  general  form: 


9“+^a-A(Y^“)+<|)=0, 

dt  dz  dz  dz 


Within  domain  (19) 


(20) 

(21) 


The  boundary  and  initial  conditions  may  be  written  as. 

y—+aa+b=0,  On  boundary  T 
dn 

o(z,0)=a0(z),  At  initial  time  zero 

In  above  equations,  o(z,t)  values  are  mean  values  of  velocity 
components,  enthalpy,  moisture  and  turbulent  correlations .♦CO 
values  are  sources/sinks  for  a.  The  Four  items  in  equa  on 
describe  the  system’s  inertia,  convection,  diffusion  and  source/sink 
strength,  respectively,  'a'  and  V  values  in  equation  20  are  used  to 

define  the  appropriate  boundary  conditions. 

A  finite-element  procedure4  was  used  for  descritization  of  spatial 
variables.  The  flow  domain  to  be  simulated  was  divided  into  small 
finite  elements.  Multiplying  both  sides  of  equation  19  by  aweighing 
factor  and  integrating  over  the  finite  element,  a  weighted  residue 
equation  set  can  be  derived.  Then,  substituting  into  the  residue 
equation  a  set  of  appropriate  interpolation  functions  for  the  weighing 
factor  and  the  element’s  variables  resulted  the  following  finite- 
element  equation  set: 


/M 


(22) 


Where  C  K_  IL  and  H,  are  the  mth-row  and  nth-column 
elements  oTsquare  matrices  which  describe  distributions  ofthe  finite 
element's  capacitance,  diffiisivity,  convectiveness  and  surface  flux, 
respectively,  and  L,  is  the  nth  element  of  a  load  vector  for  the  finite- 
element  source/sink  terms.  Finally,  assembling  the  finite-element 
equation  22  over  the  whole  flow  field  yielded  a  set  of  ordinary 
differential  equations: 


[q^I+[A:+f/+//]{A)+[i]=o 

dt 


(23) 


Where  the  vector  {A}  contains  node-points  variables  lying  within 
the  flow  domain  Q,  [L]  is  the  global  load  vector,  and  the  matrices 
rci  IK1  [U1  and  [H]  represent  global  distributions  of  capacitance, 
diffusiveness,  convection  and  surface  conductance,  respectively. 
The  finite-element  procedure  described  above  has  transformed  the 

initialboundary  valuepartial differential  equations  1  through  18  into 

«  _ _ j _ _  nrHinnrv  differential  equation  set  15. 


inmai  oouiiaary  - - .  -- 

a  large-order  system  ofthe  ordinary  differential  equation  set  23. 
Solution  of  23  may  start  from  a  Taylor  senes  expansion  of  {A} 
about  time  t„: 

ud-eS  (24) 


0tAl  , 

<AL,=<A>M0- 


dt 


dt 


where  0  is  an  implicit  coefficient.  Multiplying  both  rides  of 

equation  24  by  capacitance  [Cl  and  using  expressions  for  time 
^  -lit  _ _ o'!  dmiotinn  9.4  becomes: 


UHq^tAU-IAU 

+AtB([K+U+H]„.llA}n.i  +{Ll«-i)  (25) 

+A/(1 -e)aK+U+Hl„{A}n+lL)n) 


The  equation  set  25  may  be  solved  for  the  vector  {A}*.,  at  time 
V.  from  the  known  { A}„  values  at  time  t  by  a  Newtmfs  iterative 
process,  and  a  computer  program  has  been  written  to  faoliMe  &e 

process.  The  resultant  computer  program  can  be  used  to  simulate 

mixing  processes  in  the  marine  atmosphere.  Presented  below  is  an 


RESULTS  AND  DISCUSSION 


transient  exchange  of  moisture  in  the  sea/air  interface,  mixing  of 
moisture  in  the  marine  atmosphere,  and  formation  ofthe  cloud  layer. 
Graphs  plotted  in  figure  1  were  reported  in  a  previous  paper,  they 
were  used  as  initial  conditions  for  this  study. 


Fig.  1 :  Initial  Velocity-Components  and  Potential-Temperature 
Values 
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The  initial  conditions  shown  in  figure  1  were  obtained  on  the 
assumptions  that  the  geostrophic  wind  in  the  x-direction  was  at 
10m/s,  and  Coriolis  parameter  f  is  equal  to  1  .Ox  1  O'4  s'1,  the  potential 
temperature  was  chosen  to  be  well-mixed  at  290K,  and  the  total 
moisture  mixing  ratio  corresponding  to  100  percent  relative 
humidity  at  top  of  the  marine  planetary  boundary  layer  (MPBL). 
Just  above  the  MPBL  top  at  height  equal  to  1  km,  there  was  an 
inversion  layer  of 400-m  thick,  in  which  temperature  was  increasing 
at  0.2  K/m  and  moisture  mixing  ratio  was  decreasing  at  0.008 
g/kg.m,  respectively.  The  sea-surface  temperature  was  allowed  to 
vary  diumally  within  the  range  of 280  to  290F.  Conditions  shown 
infigure  1  were  for  the  instance  when  the  sea  surface  was  at  290  K 
For  the  simulation  run  using  initial  conditions  shown  in  figure  1 , 
it  was  assumed  that  the  sea  surface  temperature  was  raised  abruptly 
to  and  then  maintained  at  298  K,  and  the  air  total  moisture  ratio  at 
the  sea  surface  was  maintained  at  the  saturation  state.  After  a 
simulation  run  over  a  period  of  forty  hours,  a  large  body  of  data  was 
generated.  Plotted  in  figures  2  are  predicted  contours  of  temporal 
mean  physical-property  values  of  air  in  the  marine  atmosphere; 
plotted  in  figures  3  are  predicted  contours  of  Reynolds-stress  and 
turbulent-flux  values  in  the  marine  atmosphere.  Many  interesting 
characteristics  of  mixing  in  the  simulation  domain  can  be  observed 
in  these  contours.  Firstly,  rapid  interaction  at  the  air/sea  interface 
can  be  observed  during  the  initial  period  of  zero  to  600  minutes.  It 
is  followed  by  a  calmer  development  of  the  marine  planetary 
boundary  layer  for  about  600  minutes.  Then,  during  the  next  600 
minutes  transfer  of  enthalpy  and  moisture  continues,  as  can  be 
observed  from  the  predicted  contours  of  the  turbulent  thermal-  and 
moisture-flux  values  shown  in  figures  3C  and  3D.  Also,  can  be 
observed  in  Fig.  2D  is  the  formation  of  cloud  starting  at  the  200th 
minute,  rapid  deepening  of  the  cloud  layer  during  the  second  interval 
of 200  to  1200  minutes,  and  slower  growth  of  the  cloud  layer  during 
the  period  of  1200  to  1800  minutes.  Finally,  a  steady  state  is 
established  at  around  the  2400th  minute. 

Plotted  in  figure  4  are  the  predicted  steady-state  conditions  at  the 
end  of  the  simulation  period  of  the  40th  hour.  These  graphs  shown 
in  figure  4  can  be  compared  with  those  in  figure  1  for  the  initial 
conditions.  Changes  that  have  been  made  during  these  forty  hours 
can  be  observed.  Firstly,  warm  and  moist  air  at  the  sea  surface  has 
resulted  in  an  unstable  mixing  layer.  It  can  be  observed  by 
comparison  of  the  mean  horizontal  velocity  components  plotted  in 
figures  1 A  and  4  A  of  the  thickening  of  the  unstable  boundary  layer 
in  figure  4A  due  to  the  intensed  turbulent  exchange.  It  can  also  be 
observed  in  figure  4H  of  the  increased  turbulent  kinetic  energy  in 
comparison  with  the  corresponding  turbulent-energy  values  plotted 
in  figure  1.  In  addition,  variations  of  profiles  for  the  mean  wind 
velocity,  mean  potential  temperature,  mean  total-moisture  ratio, 
mean  liquid-moisture  ratios,  Reynolds-stress,  turbulent  thermal  flux, 
turbulent  moisture  flux  and  turbulent  kinetic  energy  can  be  observed 
in  figure  4. 

CONCLUSION 

In  conclusion,  a  numerical  study  of  turbulent  mixing  of  moist  air 
in  the  marine  atmosphere  has  been  made.  Results  of  predictions 
have  led  us  to  a  better  understanding  of  the  turbulent  mixing  process 
in  fire  atmosphere.  However,  complexity  of  the  turbulent  model  has 
limited  our  study  to  a  one-dimensional  model.  Physics  learned  in 
this  study  has  led  us  to  the  development  of  a  multidimensional  large- 


Tfe*. 


Fig.  2:  Predicted  Contours  of  Mean  Atmospheric  Quantities 
2A  -  Mean  Wind  Velocity,  m/s 
2B  -  Mean  Potential  Temperature,  K 
2C  -  Mean  total  moisture  mixing  ratio,  g/kg 
2D  -  Mean  liquid  moisture  mixing  ratio,  g/kg 
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Fig.  3:  Predicted  Contours  of  Turbulent  Fluxes 
3A  -  Turbulent  kinetic  energy,  mA2/sA2 
3B  -  Principal  Reynolds  stress,  mA2/sA2 
3C  -  Turbulent  thermal  flux,  m.K/s 
3D  -  Turbulent  moisture  flux,  m.kg/kg.s3 


Fig.  4:  Predicted  Mean  Quantities  and  Turbulent  Fluxes 
4A  -  Mean  velocity;  4B  -  Reynolds  stress; 

4C  -  Potential  temperature;  4D  -  Thermal  flux; 

4E  -  Total  moisture  ratio;  4F  -  Moisture  flux; 

4G  -  Liquid  moisture  ratio;  4H  -  Turbulent  energy. 
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eddy  turbulent  model  for  investigating  stability  of  the  marine  cloud 
layer.5 
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ABSTRACT 

Goals  of  our  atmospheric  research  at  University  of  the  District 
of  Columbia  (UDC)  are  to  identify  physical  processes  that 
determine  the  dynamics  of  the  cloud  layers  and  to  quantify  the  roles 
of  turbulence,  convection  and  thermal  radiation  that  play  in 
formation,  dissipation  and  stability  of  the  cloud  layers.  Our 
immediate  objectives  are  to  advance  theoretical  models,  use 
efficient  numerical  schemes  and  develop  computer  programs  to 
simulate  the  marine  cloud  layers.  Comparison  of  computer  results 
with  published  field  observations  will  yield  insights  into  the  cloud- 
layers’  physical  processes.  While  a  companion  paper  FEDSM98- 
4809  develops  a  second-order  closure  model  for  simulating 
turbulent  mixing  of  the  moist  marine  atmosphere,  this  paper  deals 
with  a  large-eddy  turbulent  model  for  simulating  stability  of  the 
marine  cloud  layers. 

INTRODUCTION 

Much  of  our  understanding  of  turbulent  processes  has  come  from 
careful  observation  and  sound  theoretical  modeling.  While  it  is  yet 
impossible  to  have  a  generalized  turbulence  theory  for  universal 
phenomena,  semi-empirical  models  with  different  degrees  of 
complexity  have  been  developed  with  confidence  in  simulating 
numerous  practical  phenomena:  the  mixing-length  theory  has  been 
used  to  model  turbulent  boundary  layers  on  flat  plates.11  the 
second-order  diffusion  model  has  been  used  for  simulating  the 
planetary  boundary  layers,7*1,9  and  large-eddy  models  have  been 
used  to  study  the  meso- scale  turbulence  in  atmospheres.  ’  This 
author  has  used  the  eddy-viscosity  model  to  simulate  a  tornado-like 
vortex,1  the  large-eddy-turbulence  model  for  end-wall  boundary 


layers  of  intense  vortices,2  the  €-k  model  for  vortex  flow  over  the 
water  surface,3,4  and  the  second-order  closure  model  for  the  marine 
cloud  layers.3,6  This  study  of  the  author  uses  a  hybrid  treatment  for 
the  atmospheric  turbulence  that  employs  a  second-order  closure 
turbulent  diffusion  model  for  simulating  atmospheric  mixing  layers 
and  a  large-eddy  turbulent  model  for  simulating  convective 
entrainment  at  the  cloud  top.  The  stratus  cloud  plays  an  important 
role  in  climatic  dynamics;  it  has  stimulated  extensive  research. 
Second-  and  higher-order  turbulent  models  have  succeeded  in 
advancing  theoretical  understanding  of  the  marine  cloud  layers. 
Complexity  of  those  models  often  makes  long-term  simulation  of 
the  cloud  layer  over  an  extensive  period  prohibitively  expensive. 
A  hybrid  model  using  second-order  bottom-up  mixing  and  large- 
eddy  top-down  convection  will  improve  in  computational 
efficiency. 

NOMENCLATURE 

Symbols: 

a,b  =  coefficients  in  boundary-layer  or  numerical  equations 

Cp  =  specific  heat 

i)  =  diffusion  parameter 

E  =  turbulent  kinetic  energy 

F  =  convection  flux  parameter 

g  =  gravitational  acceleration 

K  =  eddy  coefficient  values 

L  =  water-vapor  latent  heat  of  vaporization 

p  =  pressure 

t  =  time 

T  =  resolved  mean  static  temperature  of  atmospheric  air 


f 
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(U,  V,W)  =  resolved  mean  wind  velocities  in  (x,y,z)  directions 
(u,v,w)  =  turbulent  fluctuating  velocities  in  (x,y,z)  directions 
(x,y,z)  =  east-,  north- and  vertical-direction  coordinates 

a  =  a  generalized  variable  in  a  numerical  scheme 
P  =  buoyant  coefficient 

r  =  difiusivity  value  in  the  generalized  differential  equation 

y  =  difiusivity  value  in  the  generalized  boundary  condition 

6  =  a  length  scale  defined  by  equation  1 8 

(AxAy,Az)  =  node  points’  distance  in  (x,y,z)  directions 
€  =  a  turbulent  energy  dissipation  rate 

r\  —  a  moisture  parameter  defined  as  (L/Cp)3Q73T 
0  =  resolved  mean  potential  temperature  =  T+(gz+LQj/Cp 

0V  =  resolved  mean  virtual  potential  temperature 

=  T(l+1.609CVQ)+g^Cp 
0  =  resolved  mean  total  moisture  mixing  ratio 

=  a+o. 

0*  =  saturation  moisture  mixing  ratio  value  at  air  temperature 

A  =  a  defined  domain  under  simulation 

A  =  turbulence-length  scales 

£  =  a  moisture  parameter  defined  as  Cp T/L 

II  -  a  defined  domain  boundary 

p  =  air  density 

o  =  turbulent  Prandtl  number 

<J>  =  source  terms  in  generalized  differential  equations 

Suffices: 

app  =  for  intermediate  approximate  value  in  an  iterative  process 
e  =  for  turbulent  kinetic  energy 
f  =  for  liquid  water 
m  =  for  momentum 
o  =  for  a  reference  state 

(PJE,W3,T)  = 

for  a  variable  node  point  in  a  domain  element’s  center 
and  its  neighboring  node  pointes  at  the  east,  west,  top 
and  bottom 
(p,e,w,b,t)  = 

for  a  velocity  node  point  on  an  edge  of  the  domain  element 
and  its  neighboring  velocity  node  pointes  at  the  east,  west, 
top  and  bottom 
v  =  for  water  vapor 
0  =  for  enthalpy 

<o  =  for  moisture 

CONSERVATION  AND  TURBULENCE  EQUATIONS 

When  these  assumptions  are  made:  (1)  the  Coriolis  force  is 
negligible,  (2)  velocity  and  temperature  of  vapor  and  liquid  moisture 
are  in  equilibrium,  and  (3)  Boussinesq  approximations  are  used,  the 
conservation  equations  for  momentum,  enthalpy  and  total  moisture 
of  atmospheric  air  can  be  written  as: 


du 

dt 


=iJL(K—)*—[K  (—+—)] 

dx  "  dx  &l  "  dz  a* 


wdu  yrdU  1  dP 

dx  dz  podz 


(i) 
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w.*-[K  (§E+iK)]+2±(K  — )+Pg(©-e ) 
dt  a*  mKdz  dxn  dz  Mdz’  m  ° 

vdW  wdW  1  dP 
dx  dz  po  dz 

(2) 

( 

dt  dx  0  dx  dz  0  dz  dx  dz 

(3) 

—=JL(K  ^.)-U—-W^- 

dt  dx  w  dx  dz  “  dz  dx  dz 

1 

(4) 

dt  dx  dz  dx  dx  dz  dz  0O  dz 

+2Km(-?+KM(—+—f+2Km(—f]-e 
"  dx  mdz  dx  m  dz  1 

1 

(7) 

dx  dz 

(6) 

In  above  equations,  the  buoyancy  terms  associated  with  the  virtual 
dry  potential  temperature  have  been  defined  as  follows: 

0v=71(l  +  1.6O9Qv-O)+^ 

cp 

(7) 

TF-S  ***-'%? 

For  Clear  Air  Layers 

(8)  ‘ 

1+1.609^39  c  L  da 
dz  1+r )  dz  c  dz 

For  Cloudy  Air  Layers 

(9)  ' 

The  rate  of  dissipation  within  the  grid  volume  and  the  subgrid  eddy 
coefficient  may  be  parameterized  through  j 

€=0.l9E3n/X 

(10) 

Km=0.5ZEmX 

(11)  j 

K  2X 

k6=— (i+JR 

Ofl  Ar. 


source/sLv  **  .^uis  for  a.  Items  in  equation  19  describe  the 

(12)  variables’  time  derivative,  convection,  diffusion  and  source/sink 
strength,  respectively,  'a'  and  V  values  in  equation  20  are  used  to 
define  the  appropriate  boundary  conditions. 

A  finite-volume  difference  scheme13  was  used  for  descritization  of 

(13)  variables.  The  flow  domain  to  be  simulated  was  divided  into  small 
rectangular  elements;  shown  in  figure  1  are  examples  of  several 
such  elements.  It  can  be  seen  in  figure  1  that  variable  values  at  node 
points  of  the  elements  are  strategically  located.  The  node  point  for 

(14)  velocity  value  is  located  in  the  middle  of  the  rectangular  edge  and 
the  node  point  for  other  variables  is  located  in  the  center  of  the 
rectangular  box. 


where  A.  is  a  minimum  of  the  Blackadar’s  length  XB,  diffusion  length 
ip,  resolvable  length  scale  Xs: 


,  _  0.3526 

(15) 

B  6+3.52 

lD=0.15(^)in 

gWy 

(16) 

Xs=(AxAxA z)m 

(17) 

f ~Emzdz 


6-  0 

(18) 

f~Emdz 

.“1 

t* 

♦®W  - 

uw 

1 

- >  #0fp  — 

f  - 

->  •«! 

Ue 

k 

Fig.  1 :  An  Element  with  Neighboring  Variables  Shown 

Using  the  velocity  values  and  values  of  a  shown  in  the  figure  and 
employing  an  upwind  logic,  the  conservation  law  may  be  applied  to 
obtain  an  expression  for  the  value  of  a  at  the  node  point  P  in  terms 
of  the  a -values  at  its  neighboring  node  points: 


And  the  terms  o0,  ou,  and  ^  ,  which  are  the  turbulent  Prandtl 

numbers  of  enthalpy,  moisture  and  turbulent-energy  transports,  are  aP*P =aB*E +aifiw+aTaT+aB:tB  +b 

equal  to  0.75, 0.75  and  1 ,  respectively. 


NUMERICAL  PROCEDURE  AND  COMPUTER 
ALGORITHM 


where  coefficient  values  ‘a’  and  ‘b*  may  be  expressed  in  terms  of 
the  diffusion  and  flux  parameters  defined  as  follows: 


The  conservation  and  turbulence  transport  equations  1  to  5 
described  above  can  be  written  in  the  following  general  form: 

a-=^AZ+Mzr[-(t/.Ar),0] 

£  Ax, 

(23) 

dt  dx  dz  dx  dx  dz  dz 

within  Domain  A 

(19) 

aw=K^  +A/act(C/^z),0] 

Ar 

(25) 

The  boundary  and  initial  conditions  may  be  written  as: 

y— +aa+6=0,  On  boundary  II 

dn 

(20) 

ar=-^— +A/ax(-(t/fAx),0] 

(26) 

a(x/,f=0)=a0(x/).  At  initial  time  zero 

(21) 

r  Ax 

aB=  *  +Mzxf({/iAx),0] 

(27) 

In  above  equations,  a(x^,t)  values  are  variable  values  of  velocity 
components  U  and  W,  potential  temperature  8,  moisture  mixing 
ratio  £5,  and  turbulent  kinetic  energy  ‘E’;  <J)(x^,t)  values  are  the 

AxA z 

ap=aE+av+aT+aB+  . 

(28) 
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b=$pAxAz 


(29) 


half,  one,  five  and  ten  hours  from  the  start  of  the  simulation  run),  ( 
dissipation  of  the  cloud  layer  can  be  observed 


In  addition,  in  solving  U  and  W  values  using  equation  22,  the 
initial  pressure  values  will  have  to  be  the  estimated  values; 
consequently,  the  equation  22  will  yield  initially  approximate 
and  Ww  values.  Improved  U  and  W  values  may  be  calculated  from 
their  approximate  values  using  equations: 


U  =U  +-^L(P  -p  ) 
‘  ^  pB> 


(30) 


U 


=u 


w,app 


(31) 


WhsW>*,*-^p>-pp)  (33) 


Substituting  U  and  W  values  in  equations  30-33  into  the  continuity 
equation  6  yields  a  set  of  linear  equations  for  pressure  values  at  the 
solution  domain’s  node  points;  they  may  be  used  to  solve  for 
improved  pressure  values.  So  the  process  may  be  repeated  to  iterate 
alternatively  for  improved  values  of  velocity,  pressure  and  other 
variables  at  the  node  points. 

Using  the  linear  equations  discussed  above  for  variable  values  U, 
W,  8,  Q,  E  and  P  at  node  points,  a  computer  program  has  been 
written  to  simulate  dynamics  of  the  marine  cloud  layers. 

RESULTS  AND  DISCUSSION 

The  above-described  computer  program  may  be  used  to  study 
stability  of  the  marine  cloud  layers.  Graphs  plotted  in  figure  2  show 
characteristics  of  a  horizontally  uniform  cloud  layer  predicted  by  a 
one-dimensional  turbulent  model.6  It  may  be  noted  that  solutions  for 
this  set  of  graphs  were  obtained  by  using  the  main  stream  wind 
velocity  U  equal  to  10  m/s.  To  investigate  effects  of  cloud-top- 
warm-air  entrainment  on  stability  of  the  cloud  layer  shown  in  figure 
2,  a  stream  function  shown  in  figure  3  will  be  superimposed  to  the 
main  flow.  In  generating  the  stream  function  shown  in  figure  3,  the 
top-down  entrainment  was  assumed  to  be  at  a  maximum  rate  of  6 
cm/s  is  at  the  top-left  comer  and  the  rate  was  reduced  sinusoidally 
to  zero  at  the  top-right  comer.  In  addition,  it  is  assumed  that 
potential  temperature  of  the  entrained  air  is  at  five  degrees 
centigrade  higher  than  that  of  the  cloud  at  the  top. 

Using  the  initial  conditions  and  entrainment  rates  described 
above,  the  present  large-eddy  simulation  computer  program  has 
been  run.  Shown  in  figure  4  are  contours  of  the  initial  steady-state 
potential  temperature  values,  total  moisture  mixing  ratio  values  and 
liquid  moisture  Mixing  values,  respectively.  Predicted  dynamic 
responses  of  the  cloud  layer’s  liquid  moisture  content  to  the  warm- 
air  entrainment  at  the  top  are  shown  in  figure  5.  From  snapshots 
shown  in  this  figure  of  the  cloud  contours  at  different  times  (i.e.,  at 


CONCLUSION 

In  a  companion  paper,6  a  second-order  closure  model  was 
developed  to  simulate  sea-to-air  heat  and  moisture  transfer  and  to 
observe  theoretically  the  growth  of  the  marine  cloud  layers.  Results 
of  predictions  have  led  us  to  a  better  understanding  of  the  turbulent 
mixing  process  in  the  atmosphere.  However,  complexity  of  the  high- 
order  closure  model  had  limited  our  study  to  the  mixing  layer 
problems  in  one  dimension.  Physics  learned  in  that  study  has  led  us 
to  the  development  of  a  multidimensional  large-eddy  turbulent 
model  presented  in  this  paper.  The  extended  theory  has  enabled  us 
to  investigate  multidimensional  effects  of  convective  entrainment  at 
the  cloud  top  on  the  dynamics  of  the  marine  cloud  layers.  These 
quantitative  results  have  provided  us  at  UDC  with  theoretical  tools 
to  identify  processes  observed  physically  in  the  marine  cloud  layers. 
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Horizontal  distance  in  kilometer* 


Fig.  2:  Characteristics  of  an  Initial  Steady-State  Cloud  Layer  h<»u»*»i  dutmc  »  kaometm 

2A  -  Mean  velocity,  2B  -  Reynolds  stress; 

2C  -  Potential  temperature;  2D  -  Thermal  flux;  Fig  4;  ^  Steidy  ,ute  ^  the  Zeroth  Hour 

2E  -  Total  moisture  ratio;  2E  -  Moisture  flux; 

2G  -  Liquid  moistureratio;  2H  -  Turbulent  K.E. 
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@  Half  hour  later 


0  100  200 


Horizontal  distance  in  kilometers 

Fig.  5:  Predicted  Effects  of  Warm-Air  Enrtrainment  on  the  Cloud  Stability 
(Plotted  contours  at  different  time  are  for  liquid  water  humidity  ratio 
values  in  grams  per  kilogram) 


Copyright  ©  1998  by  ASME 


56 


Reprints  of  a  Paper  From 

1999  ASME  Fluids  Engineering  Division  Summer  Meeting,  July  18-22, 1999,  San  Francisco,  CA 
CD-ROM  Proceedings  Paper  No.  FEDSM99-7272 


APPENDIX  P 

A  HYBRID  TURBULENCE  MODEL  FOR  SIMULATING  THE  DYNAMICS 
OF  THE  MARINE  CLOUD  LAYERS 


57 


58 


Proceedings  of  the 
3Fd  ASME/JSME  Joint  Fluids  Engineering  Conference 
July  18-23, 1999,  San  Francisco,  California 


Paper  No.  FEDSM99-7272 


A  HYBRID  TURBULENCE  MODEL  FOR  SIMULATING  THE 
DYNAMICS  OF  THE  MARINE  CLOUD  LAYERS 


Joseph  Chi 

Department  of  Mechanical  Engineering 
University  of  the  District  of  Columbia 
Washington,  DC  20008-1174 
U.S.A. 


Phone:  (301)340-9794 
Fax:  (301)340-6556 
E-mail:  jchil23@aol.com 


ABSTRACT 

Goals  of  the  atmospheric  research  at  the  University  of  the 
District  of  Columbia  (UDC)  are  to  identify  physical  processes  that 
determine  the  dynamics  of  the  marine  cloud  layers  and  to  quantify 
the  roles  of  turbulence,  convection  and  thermal  radiation  that  play 
in  formation,  dissipation  and  stability  of  the  cloud  layers.  Our 
immediate  objectives  are  to  advance  theoretical  models,  use 
efficient  numerical  schemes  and  develop  computer  programs  to 
simulate  the  marine  cloud  layers.  For  these  objectives  two 
turbulence  models,  using  the  second-order-closure  and  the  large- 
eddy-simulation  (LES)  models,  respectively,  were  firstly 
formulated.  In  order  to  retain  details  of  the  second-order 
turbulence  and  calculation  efficiencies  of  the  LES  model,  a  hybrid 
turbulence  model  has  been  assembled  and  it  is  presented  in  this 
paper.  The  model  uses  a  multidimensional  framework  of  the  LES 
model  but  uses  a  set  of  unidirectional  second-order  equations  to 
calculate  the  turbulent  intensity  values.  Results  of  predicted  cloud 
profiles  using  the  hybrid  model  are  compared  with  experimental 
data  and  predictions  using  the  LES. 

INTRODUCTION 

Much  of  our  understanding  of  turbulent  processes  has  come  from 
careful  observation  and  sound  theoretical  modeling.  While  it  is  yet 
impossible  to  have  a  generalized  turbulence  theory  for  universal 
phenomena,  semi-empirical  models  with  different  degrees  of 
complexity  have  been  developed  with  confidence  in  simulating 


numerous  practical  phenomena:  the  mixing-length  theory  has  been 
used  to  model  turbulent  boundary  layers  on  flat  plates.12  the 
second-order  diffusion  model  has  been  used  for  simulating  the 
planetary  boundary  layers,8,9,10  and  large-eddy  models  have  been 
used  to  study  the  meso-scale  turbulence  in  atmospheres.11  This 
author  has  used  the  eddy-viscosity  model  to  simulate  a  tornado-like 
vortex,1  the  large-eddy-turbulence  model  for  end-wall  boundary 
layers  of  intense  vortices,2  the  €-k  model  for  vortex  flow  over  the 
water  surface,3,4  and  the  second-order  closure  model  for  the  marine 
cloud  layers.5  In  the  previous  year,  this  author  developed  a  large- 
eddy  model7  for  simulating  convective  entrainment  at  the  cloud  top 
and  a  second-order-turbulence  model6  to  study  turbulence  mixing 
processes  in  the  marine  atmosphere.  For  high  resolution,  a  detailed 
second-order-closure  model  of  turbulence  may  be  used;  for 
calculation  efficiency,  a  LES  model  may  be  used.  For  the  purpose 
of  simulating  stability  of  the  marine  cloud  layers,  it  is  necessary  to 
include  the  multidimensional  effects.  A  hybrid  turbulence  model 
that  is  built  upon  a  framework  of  the  LES  model  but  uses 
turbulence  values  predicted  by  a  second-order-closure  model  has 
now  been  developed  to  predict  stability  of  multi  dimensional 
marine  cloud  layers.  This  hybrid  model  preserves  details  of 
turbulence  but  eliminates  the  need  of  evaluating  the  multi 
dimensional  turbulence  equations,  and  It  improves  prediction 
efficiency  and  conserve  computational  resources.  Presented  in  the 
next  two  sections  of  this  paper,  respectively,  are  a  LES-based 
framework  for  the  hybrid  model  and  a  set  of  second-order- 
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turbulence  equations  for  calculating  the  hybrid  model’s  turbulence 
values.  It  is  followed  by  showing  a  numerical  scheme  that  has  been 
used  for  solving  the  model  equations.  The  predicted  cloud  profiles 
will  be  compared  with  experimental  data  retrieved  from  reliable 
web  sites  and  predictions  by  the  LES 


NOMENCLATURE 

Cp  =  specific  heat 

A,  B,  C,  D,  E,  &  F  values  in  turbulent  flux  equation 

=  empirical  constants  in  2nd-order  equations 
D  =  diffiision  coefficient  in  LES  model  equations 
f  =  geostrophic  Coriolis  coefficient 
Fr  =  thermal  radiation  flux 

g  =  gravitational  acceleration 

K  =  eddy  viscosity  and  diffusion  coefficients 
L  -  water-vapor  latent  heat  of  vaporization 

P  =  pressure 

q  =  square-root  value  of  squared  mean  of  the  turbulent 
fluctuating  velocity 


t 


=  time 


T  =  static  temperature  of  atmospheric  air 

(u,v,w)  =  turbulent  fluctuating  velocities  in  (x,y,z)  directions 
(U,V,W)  =  resolved  mean  wind  velocities  in  (x,y,z)  directions 
(x,y,z)  =  east,  north  and  vertical-direction  coordinates 

a  -  a  generalized  variable  in  a  numerical  scheme 
P  =  buoyant  coefficient 

r  =  diffusivity  value  in  the  generalized  differential  equation 

Y  ~  diffusivity  value  in  the  generalized  boundary  condition 

8  =  a  length  scale  defined  by  equation  1 8 
(Ax,Ay,Az)  =  node  points’  distance  in  (x,yyz)  directions 
e  =  a  turbulent  energy  dissipation  rate 
*1  “  a  moisture  parameter  defined  as  (L/Cp)dQ7dT 

0  —  turbulent  fluctuating  value  of  0 

0V  =  turbulent  fluctuating  value  of  0V 

0  =  mean  potential  temperature  =  T+Cgz+LClJ/C 

©V  -  mean  virtual  potential  temperature 

=  T(l+1.609Qv-Q)+gz/Cp 

k  =  Prandtl  mixing  length  constant 

\  =  turbulent  length  scale 

£  =  a  moisture  parameter  defined  as  CpT/L 
p  =  air  density 

a  =  turbulent  Prandtl  number 

0)  —  turbulent  fluctuating  value  of  Q 

£1  =  mean  total  moisture  ratio 

Q,  =  mean  liquid  moisture  ratio 

fly  =  mean  vapor  moisture  ratio 

Overlined  values 


=  turbulent-flux  values 
Other  symbols  in  numerical  procedures 
=  those  defined  in  the  text 


Suffices: 

e  =  for  turbulent  kinetic  energy 

I  =  for  liquid  water 

m  =  for  momentum 

o  =  for  a  reference  state 

v  =  for  water  vapor 

0  =  for  enthalpy 

0)  =  for  moisture 


FRAMEWORK  FOR  A  HYBRID  SIMULATION  MODEL 

When  these  assumptions  are  made:  (1)  the  Coriolis  force  is 
negligible,  (2)  velocity  and  temperature  of  vapor  and  liquid 
moisture  are  in  equilibrium,  and  (3)  Boussinesq  approximations 
are  used,  the  conservation  equations  for  momentum,  enthalpy  and 
total  moisture  of  atmospheric  air  can  be  written  as: 


EL=2±iKm^±[K 

dt  dx  mdx  dzl  m  dz  dx’1 
-U^-W—-  1  dP 


dx 


dz  p  odz 


a) 


un  U  ,vu 


—)\+2  ^(K—y^e-e ) 

dx  a?  mdz  ° 


dx  dz  o  dz 


(2) 


^=A(Ke^)+JL(Ke—)-U—-W'®  (3) 

dt  dx  dx  dz  0  dz  dx  dz 


dQ_=± 

dt  dx 


(kJQ)+JL(ku—)-u—  -  w— 

adx  dz  " dz  dx  dz 


(4) 


dU_+dW 
dx  dz 


=0 


(5) 


In  above  equations,  the  buoyancy  terms  associated  with  the  virtual 
dry  potential  temperature  have  been  defined  as  follows: 


©v=r(l+1.609Qv-fi)+i^. 


^“(0.6095-1)^2 

dz  dz  C  dz 

p 

For  Clear  Air  Layers 


(6) 


(7) 


1+1.6095^36  ,  L_dQ 
dz  1  +T|  dz  dz 

p 

For  Cloudy  Air  Layers 


(8) 


The  rate  of  dissipation  within  the  grid  volume  and  the  subgrid  eddy 
coefficient  may  be  parameterized  through 


€=0.067 q3/k 
Xn=0.41qA 
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(10) 


(22) 

dt  dz  dz  dz  X 


K  2X 

Ke=^(U^) 

°e  *>s 


(ID 


^l±Aq\—)-w—-^^--E2-uQ 
dt  dz  dz  dz  dz  X 


(23) 


K  2X 

K— (^T1) 


(12) 


^.=l(Aqk—)-^i—-w9—-E^vQ 
dt  dz  dz  dz  dz  X 


(24) 


K=- 

ea 


(13) 


(25> 

dt  dz  dz  dz  X 


where  X  is  a  minimum  of  the  Blackadar’s  length  XB,  diffusion 
length  XD,  resolvable  length  scale  Xs: 


ZLlv**) -ssS-rfe* 

dt  dz  dz  dz  A 


(26) 


,  0.35z6 

(14) 

1  6.3.52 

Wf)" 

(15) 

A.5=(AxAxAz)1/3 

(16) 

l°°qzdz 

6-Jo 

[~qdz 

Jo 

(17) 

And  the  terms  oe,  ou,  and  J ,  which  are  the  turbulent  Prandtl 
numbers  of  enthalpy,  moisture  and  turbulent-energy  transports,  are 
equal  to  0.75, 0.75  and  1,  respectively. 

SECOND-ORDER  TURBULENCE  EQUATIONS 

In  above  equations,  the  value  q  that  is  equal  to  root-mean-square 
of  the  turbulent  fluctuation  velocities  can  be  calculated  by  the 
second-order-closure  turbulent  transport  equations:9 


^=^AqX*£)-2w^-BZ{«2-£)-D£ 

dt  dz  dz  dz  X  3  3X 


dt  dz  dz  dz  X  3  3A 


dt  dz  dz  X  j  3a 

^±(2AqX^)+^<^-C^)^-Bl™ 
dt  dz  dz  dz  A 


(18) 


(19) 


(20) 


(21) 


^l+AqX^yw^-w^-Eluu 
dt  dz  dz  dz  dz  X 


(27) 


dt  dz  dz  dz  dz  X 


^=±J2AqX^yw2—^g^-E2-™> 
dt  dz  dz  dz  X 

^.=l(AqX—)-2^^-F^ 
dt  dz  dz  dz  X 


M“=±(AqX^.)-™—-^^-FldZ 

dt  dz  dz  dz  dz  X 


(30) 


(31) 


In  above  equations.  A,  is  a  characteristic  length  scale  that  is  equal 
to  the  value  of  the  Blackadar's  or  the  diffusion-length  scale  -  AB  or 
XD  -  whichever  is  the  smallest. 

In  equations  18  through  31,  die  time-derivative  terms  on  the  left- 
hand  side  model  the  transient  variation  of  turbulence  correlations. 
The  second-order  derivative  terms  on  the  right-hand  side  model 
turbulent  diffusion.  While  the  production  of  turbulence  due  to 
buoyancy  is  modeled  by  terms  with  buoyant  coefficient  p,  the 
production  of  turbulence  due  to  friction  is  modeled  by  products  of 
second-order  correlations  and  gradients  of  mean  variables.  Terms 
with  coefficients  B,  C,  E  and  F  represent  the  turbulent 
redistribution.  Turbulent  dissipation  is  modeled  by  terms  with 
coefficient  D.  Coefficients  K,  A,  B,  C,  D,  E  and  F  have  been 
determined  semi-empirically,3  having  the  values  equal  to  0.35, 
0.21, 0.46, 0.053, 0.132, 0.44,  and  0.23,  respectively. 

NUMERICAL  PROCEDURE  AND  COMPUTER 
ALGORITHM 
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The  conservation  and  turbulence  transport  equations  described 
above  can  be  written  in  the  following  general  form: 


<*> 

within  Domain  A 


The  boundary  and  initial  conditions  may  be  written  as: 


the  diffusion  and  flux  parameters  defined  as  follows: 

aE=^+Max[-(UeAz\0]  (36) 

Ax. 

aw=^+Max[{UJ^m  (37) 


Y^.+aa+fc=0,  On  boundary  II  (33) 

a(x/,r=u)=an^j.  At  initial  time  zero 

In  above  equations,  a(x^,t)  values  are  variable  values  of  velocity 
components  U  and  W,  potential  temperature  0,  moisture  mixing 
ratio  Q,  and  turbulent  stress  and  flux  values;  <J>(x^,t)  values  are 
the  source/sink  strengths  for  a.  Items  in  equation  32  describe  the 
variables’  time  derivative,  convection,  diffusion  and  source/sink 
strength,  respectively,  ’a’  and  V  values  in  equation  33  are  used  to 
define  the  appropriate  boundary  conditions. 

A  finite-volume  difference  scheme6  was  used  for  descritization 
of  variables.  The  flow  domain  to  be  simulated  was  divided  into 
small  rectangular  elements;  shown  in  figure  1  are  examples  of 
several  such  elements.  It  can  be  seen  in  figure  1  that  variable 
values  at  node  points  of  the  elements  are  strategically  located.  The 
node  point  for  velocity  value  is  located  in  the  middle  of  the 
rectangular  edge  and  the  node  point  for  other  variables  is  located 
in  the  center  of  the  rectangular  box. 


aT=^—+Max[-(UAx),0\  (38) 

A z. 


aB=^+Max[(UbAx)fi]  (39) 

A zu 


(40) 


b^tyjAxAz  (41) 


In  addition,  in  solving  U  and  W  values  using  equation  35,  the 
initial  pressure  values  will  have  to  be  the  estimated  values; 
consequently,  the  equation  35  will  yield  initially  approximate  U^p 
and  W#pp  values.  Improved  U  and  W  values  may  be  calculated 
from  their  approximate  values  using  equations: 


e 


ent  with  Neighboring  Variables  Shown 


(42) 

u«=u™+ ZiP”~Pp) 

(43) 

g- 

(44) 

A 

£1 

m 

(45) 

Using  die  velocity  values  and  values  of  a  shown  in  the  figure  and 
employing  an  upwind  logic,  the  conservation  law  may  be  applied 
to  obtain  an  expression  for  the  value  of  a  at  the  node  point  P  in 
terms  of  the  a-values  at  its  neighboring  node  points: 


aI&p-aEaE+aifaw+a7aT+aBctB+b  (35) 

where  coefficient  values  ‘a’  and  ‘b’  may  be  expressed  in  terms  of 


Substituting  U  and  W  values  in  equations  42-45  into  the  continuity 
equation  5  yields  a  set  of  linear  equations  for  pressure  values  at  the 
solution  domain’s  node  points;  they  may  be  used  to  solve  for 
improved  pressure  values.  So  the  process  may  be  repeated  to 
iterate  alternatively  for  improved  values  of  velocity,  pressure  and 
other  variables  at  the  node  points. 

Using  the  linear  equations  discussed  above  for  variable  values 
U,  W,  0,  Q,  P  and  otter  turbulence  quantities  can  be  calculated  to 
simulate  dynamics  of  the  marine  cloud  layers. 
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RESULTS  AND  DISCUSSION 

The  above-described  simulation  model  may  be  used  to  study 
stability  of  the  marine  cloud  layers.  To  increase  confidence  in  the 
model.  The  experimental  data  on  several  internet  sites  have  been 
downloaded.  Data  sources  such  as  GCE  (Goddard  Cumulus 
Ensemble),  GDACC  (Goddard  Distributive  Active  Archive 
Center),  FFI  (First  ISLCP  Field  Experiment)  and  TOGA-COARE 
(Tropical  Ocean  Global  Atmosphere  -  Coupled  Atmosphere 
Response  Experiment)  were  reviewed  to  determine  the  time 
interval,  grid  size  and  altitude  criteria.  It  was  determined  that  we 
should  secure  data  with  time  and  grid  intervals  at  six  hours  and  one 
mile,  respectively.  Two  sources  of  data  TOGA-COARE14 
(Webster  and  Lukas,  1992)  and  GCE13  (Tao,  1993)  were  then 
considered.  The  database  obtained  from  the  GCE  model  was  large. 
Eight  six-hours  data  were  organized  in  eight  different  data  files, 
each  file  containing  one  six-hour  interval  of  data.  The  data  files 
were  transmitted  via  electronic  mails. 

For  graphic  presentation,  these  data  were  decoded  at  UDC  by 
several  Fortran  programs,  and  were  translated  into  formats 
recognizable  by  data  plotting  software  packages  -  Surfer  and 
Grapher  -  which  are  available  in  the  author’s  laboratory.  Plotted 
below  are  sample  data  values  over  a  domain  with  a  height  of  3- 
km  and  a  horizontal  distance  of  500-km,  and  over  a  time  period  of 
24  hours  at  12-hours  intervals:  Fig.  2  shows  contours  of 
experimental  wind  stream  functions,  Fig.  3  sensible  and  latent  heat 
flux  values  at  the  sea/air  interface  sensible,  and  Fig.  4  the  initial 
cloud  profile  at  the  zeroth  hour. 
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Fig.  2  Contours  of  empirical  wind  stream  function 


3  Sea  surface  sensible  and  latent  heat  flux  values 


Fig.  4  Empirical  cloud  profile  at  the  zeroth  hour 


Using  the  empirical  stream  functions  shown  in  Fig.2,  the  boundary 
heat  flux  values  shown  in  Fig.3  and  the  initial  cloud  profile  shown 
in  Fig.4,  the  hybrid  model  described  in  this  paper  has  been  used  to 
simulate  the  cloud  stability.  Results  of  cloud  profiles  at  the  12th 
and  the  24th  hours,  respectively,  are  shown  in  Fig.5.  For 
comparison,  three  sets  of  data  shown  in  Fig.5  were  obtained  from 
the  different  sources:  the  experiments,  the  LES  predictions,  and  the 
present  hybrid-model  predictions.  Similarities  among  these  three 
sets  of  cloud  profiles  can  the  observed.  However,  critical 
comments  on  these  data  will  not  be  made,  until  after  the  radiation 
and  precipitation  effects  on  cloud  stability  are  reviewed  in  the 
future. 
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5A:  Experiments 


DISTANCE  IN  100KM 

SB:  LES  Model 


5C:  Hybrid  Model 


I 

Fig.  5  Comparison  of  cloud  profiles  plotted  with  data  from  three  different  sources 


CONCLUSION 

For  goals  of  research  at  UDC  to  identify  physical  processes 
that  determine  the  dynamics  of  marine  cloud  layers  and  to  quantify 
roles  of  turbulence,  convection  and  thermal  radiation  that  play  in 
formation,  dissipation  and  stability  of  the  marine  cloud  layers,  the 
author  and  his  co-workers  at  UDC  have  endeavored  to  achieve  the 
objectives: 

•  To  advance  turbulence  models  using  efficient  numerical 
schemes, 

•  To  develop  computer  simulation  programs  for  predicting 
stability  of  marine  cloud  layers, 

•  To  retrieve  web-site  data  of  marine  cloud  layers  from  reliable 
sources, 

•  To  establish  confidence  in  models  through  experimental 
verification,  and 

•  To  yield  insights  into  the  cloud’s  physical  processes. 

For  these  objectives,  several  turbulent  models,  using  the  second- 
order-closure  and  the  LES  model,  respectively,  have  been 
formulated.  A  finite-volume  procedure  has  been  developed  to 
solve  the  resultant  equations.  In  this  paper,  a  hybrid  model  has 
been  developed.  It  retains  turbulence  details  of  the  second-order- 
closure  model  and  calculation  efficiencies  of  the  LES  model. 
Results  of  predictions  with  the  hybrid  model  were  compared  with 
experimental  data  downloaded  retrieved  from  the  internet  sites  and 
with  predictions  by  LES.  Similarities  among  these  three  sets  of 
cloud  profiles  can  the  observed  in  Fig.  5.  Critical  review  on  these 
data  will  be  made,  when  effects  and  radiation  and  precipitation 
effects  on  cloud  stability  are  reviewed  in  the  future. 
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ABSTRACT 

Goals  of  the  atmospheric  research  at  die  University  of  the 
District  of  Columbia  (UDC)  are  to  identify  physical  processes 
that  determine  the  dynamics  of  the  cloud  layers  and  to  quantify 
the  roles  of  turbulence,  convection  and  thermal  radiation  that 
play  in  formation,  dissipation  and  stability  of  the  cloud  layers. 
Our  immediate  objectives  are  to  advance  theoretical  models,  use 
efficient  numerical  schemes  and  develop  computer  programs  to 
simulate  the  marine  cloud  layers.  Comparison  of  computer 
results  with  published  field  observations  will  yield  insights  into 
die  cloud’s  physical  processes.  In  this  paper  a  large-eddy  model 
for  simulating  die  marine  cloud  layers  is  presented,  experiences 
in  tracking  the  relevant  experimental  data  from  the  internet  sites 
are  described,  and  finally,  confidence  in  the  simulation  model  is 
established  by  comparing  simulation  results  with  exterimental 
data  from  reliable  sources. 

INTRODUCTION 

Much  of  our  understanding  of  turbulent  processes  has  come 
from  careful  observation  and  sound  theoretical  modeling.  While 
it  is  yet  impossible  to  have  a  generalized  turbulence  theory  for 
universal  phenomena,  semi-empirical  models  with  different 
degrees  of  complexity  have  been  developed  with  confidence  in 
simulating  numerous  practical  phenomena:  the  mixing-length 
theoty  has  been  used  to  model  turbulent  boundary  layers  on  flat 
plates."  the  second-order  diffusion  model  has  been  used  for 
simulating  the  planetary  boundary  layers, 7,8,9  and  large-eddy 
models  have  been  used  to  study  the  meso-scale  turbulence  in 
atmospheres.""7  This  author  has  used  the  eddy-viscosity  model 


to  simulate  a  tomado-like  vortex,1  the  large-eddy-turbulence 
model  for  end-wall  boundary  layers  of  intense  vortices,2  the  e-k 
model  for  vortex  flow  over  the  water  surface,3,4  and  the  second- 
order  closure  model  for  die  marine  cloud  layers.5  In  die  previous 
year,  this  author  developed  a  large-eddy  model6  for  simulating 
convective  entrainment  at  the  cloud  top.  An  effort  has  since  been 
made  by  the  author  to  extend  the  simulation  model  to  study 
stability  of  the  marine  cloud  layers.  Presented  below  are  a 
description  of  the  simulation  model,  a  collection  of  experimental 
data  from  reliable  internet  sites,  and  a  comparison  of  computer 
simulation  results  with  experiments  to  establish  confidence  in  die 
simulation  model. 

NOMENCLATURE 

Symbols: 

Cp  =  specific  heat 

D  =  diffusion  parameter 

E  =  turbulent  kinetic  energy 

F  =  convection  flux  parameter 

g  =  gravitational  acceleration 

K  =  eddy  coefficient  values 

L  =  water-vapor  latent  heat  of  vaporization 

P  =  pressure 

t  =  time 

T  =  resolved  mean  static  temperature  of  atmospheric  air 
(U,V,W)  =  resolved  mean  wind  velocities  in  (x,y,z)  directions 
(x,y,z)  =  east-,  north- and  vertical-direction  coordinates 

P  =  buoyant  coefficient 

r  =  diffUsivity  value  in  die  generalized  differential 
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equation 


i 


(x,y,z)  =  east-,  north-  and  vertical-direction  coordinates 

P  =  buoyant  coefficient 

T  =  diffiisivity  value  in  the  generalized  differential  equation 

y  =  diffiisivity  value  in  the  generalized  boundary  condition 

6  =  a  length  scale  defined  by  equation  1 8 

(Ax,Ay,Az)  =  node  points’  distance  in  (x,y,z)  directions 
e  =  a  turbulent  energy  dissipation  rate 

T)  =  a  moisture  parameter  defined  as  (L/CP)3Q73T 
9  =  resolved  mean  potential  temperature  =  T+fgz+LQ,  )/Cp 

8„  =  resolved  mean  virtual  potential  temperature 
«  T(l+1.609Qv-Cl)+gz/Cp 
Q  =  resolved  mean  total  moisture  mixing  ratio 
=  0,  +  Q, 

Q‘  =  saturation  moisture  mixing  ratio  value  at  air  temperature 

k  =  turbulence-length  scales 

E,  =  a  moisture  parameter  defined  as  CPT/L 

p  =  air  density 

a  =  turbulent  Prandtl  number 

Suffices: 

e  =  for  turbulent  kinetic  energy 

(  =  for  liquid  water 

m  =  for  momentum 

o  =  for  a  reference  state 

v  =  for  water  vapor 

0  =  for  enthalpy 

to  =  for  moisture 

CONSERVATION  AND  TURBULENCE  EQUATIONS 

When  these  assumptions  are  made:  (1)  the  Coriolis  force  is 
negligible,  (2)  velocity  and  temperature  of  vapor  and  liquid 
moisture  are  in  equilibrium,  and  (3)  Boussinesq  approximations 
are  used,  die  conservation  equations  for  momentum,  enthalpy  and 
total  moisture  of  atmospheric  air  can  be  written  as: 


—=2-^K—)+—[K(—+^L)] 

dt  dx  mdx  dz  "  dz  dx 


L,dU_wdU_  1  dP 
dx  dz  p  dz 


(i) 


+4-Ike^)  +[-| 


ae 


dt  ‘  dx  dzJ  dx ‘  Edx  dz  Edzj  L  eo  a 
”  dx  "  dz  dx  mdz 


(5) 


I 


cbc  dz 


(6) 


In  above  equations,  die  buoyancy  terms  associated  with  the  virtual 
dry  potential  temperature  have  been  defined  as  follows: 

9V=7(1+1.609QV-Q)+^-  (7) 


££l=^®+(0.609?-l)— — 
dz  dz  Cp  dz 

For  Clear  Air  Layers 


(8) 


i 


d8y_  l+l-609gT)d9_^  L  dQ 

dz  1  +T|  dz  Cp  dz 
For  Cloudy  Air  Layers 


(9) 
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The  rate  of  dissipation  within  the  grid  volume  and  die  subgrid 
eddy  coefficient  may  be  parameterized  through 

e=0.19 Eia/k  (10) 


i 


K=0.5SEiak  01) 

fit 

\ 


2A(^)+pg(e-eo) 

dt  dx1  m  dz  dxn  dz  mdz;  ™  °* 

_| jjdW  wdW_  1  dP 
dx  dz  p  dz 
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ao_a  ae,  a  ae,  rrae  TJ/ae 

dt  dx  0  dx  dz  6  dz  dx  dz 
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dt  dx  adx  dz  " dz  dx  dz 
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where  A,  is  a  minimum  of  the  Blackadar’s  length  XB,  diffusion 
length  XD,  resolvable  length  scale  A.s: 


,  0.35z6 

'  6*3.51 

(15) 

g°  °v 

(16) 

As=(AxAxAz)l/3 

(17) 

("Emzdz 
fi-  0 

r Em<t 

(18) 
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And  die  terms  o9,  au,  and  &  ,  which  are  the  turbulent  Prandtl 
numbers  of  enthalpy,  moisture  and  turbulent-energy  transports,  are 
equal  to  0.75, 0.75  and  1,  respectively. 

INTERNET  ATMOSPHERIC  DATA  TRACKING 

Several  data  sources  such  as  Goddard  Cumulus  Ensemble 
(GCE),  GDACC  (Goddard  Distributive  Active  Archive  Center), 
FIFI  (First  ISLCP  Field  Experiment)  and  TOGA-COARE 
(Tropical  Ocean  Global  Atmosphere  -  Coupled  Atmosphere 
Response  Experiment)  were  reviewed  to  determine  the  time 
interval,  grid  size  and  altitude  criteria.  It  was  determined  that  we 
should  secure  data  with  time  and  grid  intervals  at  six  hours  and 
one  mile,  respectively.  Two  sources  of  data  GCE13  and  TOGA- 
COARE14  were  then  considered. 

The  scientific  goals  of  COARE  are  to  describe  and  understand: 
(1)  the  principal  processes  responsible  for  the  coupling  of  the 
ocean  and  the  atmosphere  in  the  western  Pacific  warm-pool 
systems,  (2)  the  principal  atmospheric  processes  that  organize 
convection  in  the  warm-pool  region,  (3)  the  oceanic  response  to 
combined  buoyancy  and  wind-stress  forcing  in  the  western  Pacific 
warm-  pool  region,  and  (4)  (he  multiple-scale  interactions  that 
extend  the  oceanic  and  atmospheric  influence  of  the  western 
Pacific  warm-pool  system  to  other  regions  and  vice  versa.  To 
carry  out  the  goals  of  TOGA-COARE,  three  components  of  a 
major  field  experiment  have  been  defined:  interface,  atmospheric, 
and  oceanographic.  The  experimental  design  calls  for  a  complex 
set  of  oceanographic  and  meteorological  observations  from  a 
variety  of  platforms  that  carry  out  remote  and  in  situ 
measurements.  The  resulting  high-quality  dataset  is  required  for 
the  calculation  of  the  interfacial  fluxes  of  heat,  momentum  and 
moisture,  and  to  provide  ground  thruth  for  a  wide  range  of 
remotely  sensed  variables  for  the  calibration  of  satellite-derived 
algorithms.  The  ultimate  objective  of  the  COARE  dataset  is  to 
improve  air-sea  interaction  and  boundary-layer  parameterizations 
in  models  of  the  ocean  and  the  atmosphere,  and  to  validate 
coupled  models. 

Internet  web-sites  data  were  used  to  review  the  TOGA-COARE 
datasets.  The  web  site,  http://kiwi.atmos.colostate.edu/scm/toga- 


coare.html  was  found  to  contain  mean  data  over  the  TOGA- 
COARE  IFA  region,  which  were  thought  appropriate  for  testing 
model.  Using  FTP  commands  the  data  files  were  transferred 
compressed  files  were  uncompressed  and  converted  to  ASCII  files 
and  transmitted  to  a  UDC  workstation  via  electronic  mails. 
However,  owing  to  their  coarse  grid  sizes,  higher  resolution 
datasets  are  required  for  the  model  refinement. 

The  Goddard  Cumulus  Ensemble  Model  (GCE)  is  maintained 
by  Mesoscale  Atmospheric  Processes  Branch  (MAPB)  at  Goddard 
Space  Flight  Center  (NASA/GSFC).  Scientists  in  NASA/GSFC 
continuously  maintain  and  upgrade  the  system,  and  the  GCE 
model  is  simulated  and  operated  by  GSFS  system  analysts  and 
operators.13  The  GEC  model  may  be  used  to  interpolate 
experimental  data  with  great  degrees  of  accuracy.  The  following 
procedures  were  followed  to  obtain  the  data  from  both  sources: 

Super  computer  CRAY  at  NASA  was  used  for  obtaining, 
organizing,  processing  the  data.  The  data  files  obtained  in  this 
process  were  transmitted  to  the  UDC  workstation  via  electronic 
mail.  For  GCE  Model  data,  scientists  at  GSFC  provided  the 
author  and  his  co-workers  at  UDC  with  the  database  and  a  Fortran 
program  including  a  list  of  parameters  via  electronic  mail.  Fortran 
programs  were  further  developed  to  modify  the  data  into 
appropriate  formats  for  the  files.  The  data  from  GCE  model  were 
processed  using  the  Fortran  programs  in  UNK  environment. 
Then,  they  were  transmitted  to  the  UDC  workstation  in  two  stages. 

The  first-stage  data  transmitted  included  the  following 


parameters: 

parameter  (nx=512,nz=22) 

parameter  (nhr=48)  !  number  of  hours  model 

simulations  (h) 

parameter  (interval=21600)  !  Interval  to  read  in  data  (s) 

parameter  (ntimes=nhr*3600/interval) 

real  latent  (nx,ntimes)  !  surface  latent  heat  flux 


real  pO  (nz) 
real  qcg(nx,nz,ntimes) 
real  qci(nx,nz,ntimes) 
real  qcl(nx,nz,ntimes) 

real  qcs(nx,nz,ntimes) 
real  qm(nx,nz,ntimes) 


(W/m**2) 

!  air  pressure  (mb) 

!  graupel  mixing  ratio  (g/kg) 

!  cloud  ice  mixing  ratio  (g/kg) 
!  cloud  water  mixing  ratio 
(g/kg) 

!  snow  mixing  ratio  (g/kg) 

!  rain  water  mixing  ratio 
(g/kg)  . 


real  rad_lw(nx,nz,ntimes) 

real  rad_sw(nx,nz,ntimes) 

real  sensible(nx,ntimes) 

real  temp(nx,nz,ntimes) 
real  u(nx,nz,ntimes) 
real  v(nx,nzjitimes) 
real  vap(nx,nz,ntimes) 

real  w(nx,nz,ntimes) 
real  zl(nz) 

real  z2(nz) 


!  long  wave  radiative  heating 
rate  (K/hr) 

!  short  wave  radiative  heating 
rate  (K/hr) 

!  surface  sensible  heat  flux 
(W/m**2) 

!  air  temperature  (K) 

!  u-wind  speed  (m/s) 

!  v-wind  speed  (m/s) 

!  water  vapor  mixing  ratio 
(g/kg) 

!  w-wind  speed  (m/s) 

!  model  grid  height  at  zl  levels 
-  for  most  variables  here  (m) 
!  model  grid  height  at  z2 
levels  -  where  2  is  (m) 
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nx=512  (512  miles,  data  grid  =  one  mile) 

nz=22  ->  pressure  levels 

ntimes=32  (6-hrs  intervals) 

The  second-state  data  transmitted  included  the  following 
parameter: 

real  tke(nx,nz)  Iturbulent  kinetic  energy 

(cm*cm/s/s) 

real  du(nx,nz)  !  u-momentum  turbulence  rate 

(cm/s/s) 

real  dv(nx,nz)  !  v-momentum  turbulence  rate 

(cm/s/s) 

real  dw(nx,nz)  !  w-momentum  turbulence  rate 

(cm/s/s) 

nz=22  ->  pressure  levels 

ntimes=32  (6-hrs  intervals) 

The  database  obtained  from  the  GCE  model  was  large.  Eight 
six-hours  data  were  organized  in  eight  different  data  files,  each 
file  containing  one  six-hour  interval  of  data.  The  data  files  were 
transmitted  via  electronic  mails. 

DATA  PRESENTATION  AND 
COMPARISON  WITH  PREDICTION  RESULTS 

For  graphic  presentation,  these  data  were  further  decoded  at 
UDC  by  several  Fortran  programs,  and  were  translated  into 
formats  recognizable  by  data  plotting  software  packages  -  Surfer 
and  Grapher  -  which  are  available  in  the  author’s  laboratory. 
Presented  below  are  plotted  sample  data  values  over  a  domain 


Fig.  1  Contours  of  empirical  wind  stream  function 


•  Contours  of  empirical  wind  flow  stream  functions  values 

•  Contours  of  empirical  total  moisture  mixing  ratio  values 

•  Contours  of  empirical  liquid  moisture  muting  ratio  values 

•  Sea/air  interface  sensible  and  latent  heat  flux  values 


Fig.  2  Contours  of  expirical  total  moisture  mixing  ratio 


Fig.  3  Contours  of  empirical  cloud  profiles 


with  a  height  of  3-km  and  a  horizontal  distance  of  500-km,  and 
over  a  time  period  of  24  hours  at  12-hours  intervals: 
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Fig.  4  Sea  surface  sensible  and  latent  heat  flux  values 


These  empirical  data  can  now  be  used  to  verify  the  simulation 
results  predicted  by  the  model  described  in  this  paper.  For  these 
purpose,  using  the  emiprical  boundary  streamfunctions  shown  in 
Fig.  1  and  die  sea-surface  latent  and  sensible  heat  fluxes  shown  in 
Fig.  4  as  boundary  conditions,  and  die  empirical  data  at  the  zeroth 
hour  as  initial  conditions,  computer  simulation  of  the  cloud 
stability  has  been  made.  Several  snapshots  of  the  predicted  cloud 
files  from  the  zeroth  to  the  24th  hour  are  plotted  in  Fig.  5  for 
comparison  with  the  empirical  cloud  profiles  shown  in  Fig.  3. 
Although  it  is  not  expected  to  have  exact  agreement  between  the 


and  precipatation  effects  in  die  theoiy,  similarity  of  transient  cloud 
profiles  in  these  two  figures  has  increased  confidence  in  the 
theoretical  model. 


CONCLUSION 

For  goals  of  research  at  UDC  to  identify  physical  processes 
that  determine  the  dynamics  of  marine  cloud  layers  and  to 
quantify  roles  of  turbulence,  convection  and  thermal  radiation  that 
play  in  formation,  dissipation  and  stability  of  the  marine  cloud 
layers,  die  author  and  his  co-workers  at  UDC  have  endeavored  to 
achieve  the  objectives: 

•  To  advance  turbulence  models  using  efficient  numerical 
schemes. 


!  1 


I 


©OTH  HOUR 


distance  in  iookm 


Fig.  5  Snapshots  of  predicted  cloud  profiles 


•  To  develop  computer  simulation  programs  for  predicting 
stability  of  marine  cloud  layers, 

•  To  retrieve  web-site  data  of  marine  cloud  layers  from  reliable 
sources, 

•  To  establish  confidence  in  models  through  experimental 
verification,  and 

•  To  yield  insights  into  the  cloud’s  physical  processes. 

For  these  objectives,  several  turbulent  models,  using  the  second- 
order-closure  and  the  LES  model,  respectively,  have  been 
formulated.  A  finite-volume  procedure  has  been  developed  to 
solve  the  resultant  equations. 

In  this  paper,  for  experimental  verification  of  the  LES  model. 


simulation  results  have  been  compared  with  relevant  empirical 
data  retrieved  from  reliable  sources.  Several  data  sources  such  as 
TOGA-COARE  (Tropical  Ocean  Global  Atmosphere  -  Coupled 
Ocean  Atmosphere  Response  Experiment),  GCE  (Goddard 
Cumulus  Ensemble),  FIFE  (The  First  ISLCP  Field  Experiment), 
GDAAC  (Goddard  DAAC)  and  others  were  evaluated  based  on 
the  time-interval,  grid-size  and  domain  criteria.  According  to 
these  criteria,  two  sources  of  data  TOGA-COARE  (experimental) 
and  GCE  model  (semi-empirical)  were  examined  in  detail. 
Through  cooperation  with  the  NASA  scientists  and  a  permission 
of  using  their  CRAY  computer,  several  sets  of  the  GCE  data  were 
decoded  and  transmitted  to  the  author.  These  data  have  been  used 
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to  verify  the  UDC’s  model  simulation  results.  Figure  5  shows 
snapshots  of  cloud  profiles  that  have  been  predicted  by  the  LES 
computer  program  under  the  same  wind-stream-function,  air/sea 
interface  sensible  and  latent  heating  fluxes  and  other  boundary- 
condition  values  that  were  provided  by  the  empirical  data. 
Although  predictions  are  not  expected  to  have  a  complete 
agreement  with  the  experimental  data 

(as  precipitation  is  not  considered  in  the  present  theory),  a  good 
degree  of  similarity  can  be  observed  between  the  predicted  cloud 
profiles  shown  in  figure  5  and  those  empirical  cloud  profiles 
shown  in  figure  3. 
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